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GENERAL METHOD OF EXACT SOLUTION OF THE 
CONCENTRATION-DEPENDENT DIFFUSION EQUATION 


By J. R. PHm1p* 
[Manuscript received November 11, 1959] 


Summary 


Only three forms of D(6) have previously been known to yield exact solutions 


of the equation 
ob a 06 
=— =~ | D(0)—}. 
ot Ox Ox 
subject to the conditions 8#=0, z>0,t=0; 0=1,x=0,t>0. The present paper reports a 


general method of establishing a very large class of D(0) functions which yield exact 
solutions. A similar method enables exact solutions of the same equation subject to 


1 
the conditions 8=0, x>0, and §6=1, x<0, t=0; | xd§=0, t>0. In this case also 
Ones 
a very large class of D(8) functions yield exact solutions. Examples are given for both ~ 


cases. 

Many of the exact solutions which are most readily found tend to lead to zero or 
infinite values of D at one or two points of the 6-range. Means of avoiding this difficulty 
are devised. Practical use of the method is discussed. : 


I. INTRODUCTION 
Hitherto, the mathematics of concentration-dependent diffusion has depended 
almost exclusively on numerical analysis. According to Crank (1956, p. 166), 
the only known “formal solutions ” in concentration-dependent diffusion are 
due to Fujita (1952a, 1952b,.1954). They provide solutions of the equation 


00. a .) 
>= seme | iy aks she a chan ce: Ne es 1.1 
51 aa Paz) (1.1) 
subject to the conditions 
0=0, w>0, t=0; ) 
0=1, #=0, t>0; J 


and are for the following D(6) functions : 
D=D,/(1—r0); D=D {A=}? ; D=D,/(1+2a0 +50?). 


Philip (1955) showed how (1.1) subject to (1.2) could be solved analytically 


for D, an n-step function of 0. However, the method is cumbersome when n 


is large. i 
- This paper presents a very general method of obtaining exact solutions of 
(1.1) subject to (1.2). It embraces all possible exact solutions, including as 


particular cases those of Fujita and those implicit in Philip (1955). The method 
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is also applicable to all possible exact solutions of (1.1) subject to the following 
conditions, which are also of practical importance : 


1 
@=0, «>0 and 0=1, «<0, t=0; i} wa0=0; 10. eee Tals) 
0 
Extension of the results of this paper to the slightly more general cases of 
(1.2) and (1.3) with 6=6,, 0=0,, in place of 6=0, 0=1, is trivial, involving only 
the introduction of the linear transformation 6,—(6—6,)/(6,—9,) and solution 
in terms of 0,. 


II. Exact SOLUTIONS AND Exact FUNCTIONS 

The use of the adjective ‘‘ exact’ to distinguish solutions obtainable in 
terms of the functions of analysis from those obtainable only numerically is well 
established (e.g. Carslaw and Jaeger 1959, p. 91). This usage is convenient and 
unambiguous, though no definition of the term seems to have been given in the 
literature. 

For the purposes of the present work, it is desirable to define an exact function 
of a real variable. We say that y is an exact function of the real variable x in 
some interval of x, provided y can be specified without approximation for all 
values of x in the given interval by means of a finite number of explicit formulae. 

Further, we may then define an exact solution as a solution expressible in 
the form of an exact function. 


III. GENERAL PRINCIPLES OF THE METHOD 
It is well known (cf. Philip 1955) that the substitution 


I oe a 3 so atta eR ee ee (3.1) 
enables (1.1) subject to either (1.2) or (1.3) to be reduced to 


0 
| od0 == 2DdG/dg, i cess. eee (3.2) 
0 
that is, 


0 
D= —idg/a0. | O10. ete. cee (3.3) 
0 


It follows that the solution of (1.1), subject to either (1.2) or (1.3), ¢(6), 
exists in exact form, so long as D(8) is of the form 


6 
D=—14F/d6. { Fag.) ee arene ae (3.4) 
0 


where F is any single-valued exact function of 6 in the interval 0<0 <1, which 
satisfies certain, not very restrictive, conditions, which we indicate in the following 
paragraphs.} 


t Note that D is expressible in terms of known functions if F (as well as satisfying these 
conditions) is integrable in terms of known functions. 
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i) Conditions imposed by Conditions governing (1.1).—Evidently, the exact 
solution for D(0) of form (3.4) is simply 9>=F. F must therefore satisfy the 
conditions to be imposed on (p in) equation (3.2). It follows that, for (1.1) 
subject to (1.2) 


? 


EEC L Joes) ees ate es ee te Ae (3.5) 
and that, for (1.1) subject to (1.3), 
VCS OE Cais ae Roe a ares (3.6) 


(ii) Conditions imposed by Requirement that D Ewxists—For D to exist for 
‘ 0 

all values of 6 in the range 0<6<1, it is necessary that | Fd0 and df/d6 exist 
0 


throughout this 6-range. (However, if a finite number of discontinuities in D 
are allowable, or if D is permitted to be infinite at a finite number of points in 
the range 0<6<1, df/d6 may either not exist or be infinite at the appropriate 
finite number of points in the 6-range.) 


(ili) Condition imposed by Requirement that D>0.—The flux of the diffusing 
entity in the sense « positive, Q, is equal to —D00/dv=—t-iDd0/do. Now 
both (1.2) and (1.3) describe phenomena in which for D>0, Q>0. It follows 
that d6/d¢ <0 in both cases. Accordingly we have, for the two cases we consider, 
the further condition 

dF jd0= 0: (0—.0 <2 1). 4 eee ae cok (3.7) 


IV. SoLuTions oF (1.1), (1.2) 
It is seen that exact solutions of (1.1) subject to (1.2) may be established 
at will, simply by selecting / functions satisfying the conditions set out in 
Section III. Some typical results of this elementary process are given in Table 1. 


Tt will be observed that the D functions in this table tend to have lim D(0) 


6-0 
either zero or infinite.* It would be a serious limitation on the generality of 
the present methods if none of our exact solutions were for D(0) non-zero, but 
finite.t Accordingly, we now give some attention to this question. 


A well-known result in linear diffusion (Carslaw and Jaeger 1959, p. 59) 
for phenomena governed by (1.2) is equivalent, from the present viewpoint, 
to the statement that 


for F—2D} inverfc 0, D=constant=D,. ...... (4.1) 


The notation inverfe denotes the inverse of erfe (Philip 1955). 


* From this point on we shall write D(0) for lim D(8). Note that we use the form 
6-0 
lim D(6) since, strictly, the value of D at 6=0 is irrelevant to the phenomena. We are con- 


6-0 : ; 
cerned only with values of D(@) in the interval 0<0<1. 

+ This is not to imply that no practical interest attaches to cases with D(0)=0. For a 
discussion of such cases and certain applications see Philip (1956, 1957). 
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We therefore put 


F=2D) inverfe §+/(6).; Do finite, non-zero me Bed 
in (3.4), and obtain 
6 
P= D, +404 (7 [ ja —B as) 4a. | fA9 pal oie eeee (4.3) 
0 
where 
| B(9)=2n-! exp [—(inverfe 6)?]. ...---- +++ seer eee (4.4) 
TABLE 1 


SOME SIMPLE CASES OF EXACT SOLUTIONS OF (1.1), (1.2) 


| 


No. D ) | Remarks 
1 4n0"{1—6"/(1+-n)] 1—6” n>0 
n 
3 —=6\e— 1/16)? = 
Toran de 6) —(1—9)2"] (1—6)” n>0 
9-n 
3 mol =] 6-"—1 0<n<l 
1—n 
4 4 sin? 4n0 cos $70 
5 4.cos 4nO[cos 4n6+ $r9—1] 1—sin $n0 
6 cos-1 8+1 
6 lero =) cos—! 6 
4n—(1—8) sin-? (10) 
u 1 —1,(1— 
8 sin-! 6 
# (0+ ya—) v(1—6") 
9 4 sin 70(70+sin 70) cos? $70 
10 +; sin? 70(5-+cos 70) cos® 470 
11 4079 [1—e”® + ne”) er—end n>0 
12 Jo —79[1—e —”9 ne —”] e—n_»—n n>0 
13 4(1—log 6) | ayia WAG. § 
a = (—log 6)"—} ieee 27 n>0. T(x,p) 
20 ' denotes the incom- 
x {I'(n+1)—T'(—log 9, n+1)} plete I'-function 
['s —%¢P —1dz,, 
0 


This result depends on the following relationships, which are readily established : 
Ms (inverfe ES — = mat inverfe 0 as me Set | 


. The properties of inverfc 0 and B(0) are treated in more dabei in a Philip (1960). 
Tables of these two functions are also given there. . . 
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Reverting to (4.3), we are concerned with the conditions under which .D(0) 
is finite, non-zero. Evidently, it is a sufficient condition that both 


Bm (if 20} ond aim (25) 


be finite (including zero). The slightly stronger condition that these limits 
be zero ensures that D(0) is equal to Do: See the imposition on-f of the 


conditions 
6 . 
iim Fal jad) ~0 pane Mts pAze (4.6) 
— s 0 
d 
_jim (75) E01 wae ere ee (4.7) 


ensures that the inverfc term in (4.2) accounts completely for the finite but 
non-zero value of D(0) in (4.3). 
It is simply shown that (4.6) and (4.7) are equivalent to the single condition ~~ 


that . 
jim eS EE MANE Pigs (4.8) 
Note that since (Philip 1960) 
ay inh a ee (4.9) 
(4.8) is equivalent to 
Le 18 1h et (4.10) 


lopli 


f is thus free to become infinite at 6=0, but it must approach infinity more 
slowly than (—log 6)!. Note that it is necessary, but not sufficient, that foo 
more slowly than §-*, where « is any non-zero positive quantity. 

Accordingly, there exists an exact solution of (1.1) subject to (1.2) with 
D(0)=D, (finite, non-zero) corresponding to every D(@) of the form (4.3) for 


which © 
(i) (2D8 inverfc +f) satisfies the appropriate conditions on F set out in 


Section II; 
(ii) f satisfies condition (4.8). 
A further result is that, for | 
F=2D} inverfc 6+f (D, finite, non-zero), wes (4.11) 


D=D,+3D)( (5 f0—B55)- i ra0-4 AF ae sup (4.12) 


For the results (4.2), (4.3), the restriction on df/d@ was simply 
dffdGE2DyBi eile Tes 46° (4,13) 
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We here require the somewhat stronger restriction 
| (afd |2oDye. eee eee (4.14) 


We note in passing that condition (3.5) reduces to the condition on f that 


It will be observed that a very large class of functions f satisfy these various 
requirements, and that many of them may be found by elementary means. Some 
typical results are presented in Table 2. It is, of course, a simple matter to 
extend this table to a much greater length; but there seems little purpose in 
enumerating further examples. 


TABLE 2 
SOME EXACT SOLUTIONS OF (1.1), (1.2) witH D(0) FINITE NON-ZERO 


| 
No. D © Remarks 


| 
67 L 
1 | Dot inor(1— eo 2D6 inverfe 6 +(1—6”) | n>0O. For minus sign 
wd | D,>7n*/16 
3 20 6” = | 
Do | =) 1— Bnor-1 
£209 | 4 l+n eek 
n 
Do+—— [(1— 6)” —1—(1—6)2"] 2D} inverfe 6 +(1—6)” | n>0. For minus sign 
2(n+1) gel 
D, has lower limit 
} | 
D 11 —6)\r +2 
+4 | ae (La O)0 =] 
+Bn(1—6)” | 
é rs 8 
3 | Do+# sin? 4n0 +4D3 sin $0: EB +7] 2Dp3 inverfe O+cos 4x6 | For minus sign 
= ; Dy> x16 
4 Dy +4e”[1—e”® + n6e”] 2D6 inverfe 8 +(e"™—e"®) | n>O. For minus sign 


Dy> mn®e2"/16 


3| 2 
+4D5 aa 1—e”® + nOe”) + Bner® 


V. SOLUTIONS oF (1.1), (1.3) 

Similarly, we may establish exact solutions of (1.1) subject to (1.3) at will, 
simply by selecting # functions satisfying the appropriate conditions stated in 
Section III. Typical results are given in Table 3. It will be noted that the 
entries 1 and 11 of the table are for D(@) functions symmetrical about 6=4. 
It is evident that these results follow immediately from Nos. 2 and 4 of Table 1, 
and that many more ‘‘ symmetrical ” exact solutions may be readily constructed 
from the results of the preceding section. 

We here again encounter the difficulty which arose previously, namely, 
that F functions chosen without certain precautions tend to lead to zero or 
infinite values of D. In Section IV, this occurred only for D(0); but here the 
trouble arises also for D(1), which we write for lim D(6). 

O1 
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In this case we may avoid the difficulty by using a known result for diffusion 
in composite media governed by conditions (1.3) (Carslaw and Jaeger 1959, 
p. 87). We find that, from the present viewpoint, this result is equivalent to 
the following statement. 


Let F be specified by (5.1) below: 


1>0>6’, F=—2D} inverfe {(1 —6)/(1—6’)} ; | 


0<0<6’, F=2Dj inverfe (6/6’) ; .. (5.1) 
where | 
6’ =D} /(D§ + D)). 


TABLE 3 


SOME SIMPLE CASES OF EXACT SOLUTIONS OF (1.1), (1.3) 


No. D co) Remarks 
1 | —"— (1 —26)"—-11—26)"] (—20)" Uhh thy, <- 
2(n+1) 
nor 1 ae ey, 
2 1—6” ode es 
2np1): n+1 
1 
a ee 0<n<l 
3 9-"—1 — 
—n)\ ) 1 : 
n(1—6)” ees PEA 
4 1—(1—0)” 
ao an 
5 n(1—6) * (18) -"—1] —_ —(1—8)~” 0<n<l 
2(1—n) —m 
4 sin $7r0(sin 4r6—8) cos ae 
4. cos 4n6[cos 4n86+0—1] 2/m—sin 476 
§ 1+6 cos-16—0 ; tren) 
Vv (1—8?) 
ey ‘au Rad aaa sin-(1—0)—jn+1 
/ (20—6?) 
10 10 a aa V(1—62)—4n 
/(1—6?) 
ll 4 sin? 70 cos 770 
i 3 tn0—4/3 
ae aul eat) cos r9—40 cosec $70] cos? $7 [37 
- n__] e?—1—nen® n may be either 
3 ts a ar : positive or 
negative 
14 —tHlog 8 —(log 9+1) 
15 ~ (log 6)"—1 (log 0)"—T'(n +1) n>0 
20° 


x {Tn +1) -(1—0) I'(—log 6, n+ 1)} 


8 Jo keene e 
Then the corresponding D function is : 


1>6>6’, D=constant=D, ; ) 


Rees (5.2) 
0<60<0’, D=constant=D,. J 
We therefore put F into the form (D,, D, finite, non-zero) 
1>0>0', F=—2D} inverfe {(1— Pe 6')} +f (6 5) (5.3) 


0<0<0’, F=2D% inverfe (6/6’) + (6 


Using (3.4), we obtain the corresponding D(@) for diffusion subject to conditions 
(1.3): 


2 af 
Sy, pt | Tene | (a0 aie 6yB(— a 


ie 
iy | a0: 
' t r , 
0<6<60', D=D,+4D5|2/0 aI “fao— —6’'B (6/6') 56 
Of, 
ae iL ‘fa0, 
sje: oon Ae ea eemla eae (5.4) 
The condition that D(0)—D, reduces to 
; f = 
pea inverfc (0/0)? og © Fae wes ate 6 (5.5) 


entirely analogous to (4.8) in the previous Section. We have similarly that the 
condition that D(1)=—D, is that 


f 


ll eer (pases (5.6) 


It follows that there exists an exact solution of (1.1) subject to (1.3) with 
D(0)=D,, D(1)=D, (both finite, non-zero) corresponding to every D(6) of the 
form (5.4) for which 

(i) F of the form (5.3) satisfies the appropriate conditions set out in 
ee Duis 
i) af Raticnas conditions (5.5) and (5.6). 


. will be noted that, in this case, when df/d6 is continuous at 6=60’, 
discontinuity of magnitude D,—.D, occurs in D(@) at this point. The Pet 
that dF/d6, and therefore D(@), is continuous at 6=0' is that 


Wire ek 


That is, a discontinuity in df/d@, defined in (5.7), is necessary if a discontinuity 
in D(8) is to be avoided at 0=6'.. It is seen that usually f will be required to 
satisfy this further condition. 


a 
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ee § RLS ale helpful to replace f by g +f’, where the discontinuity in df/dé 
at §=@’ is contained in dg/d6, with df’/d6 continuous at this point. Condition 
(3.6) on F reduces here to the requirement that 


jal 
| fa0—0, 
0 


so that it is convenient to take g such that 


The family of linear g functions satisfying (5.8) is 


1>6>6', g=a+b(0—0’); ) 
0<0<0', g=a+(b—c)(6—6’). | 


Here a is an arbitrary constant, 


b=(Di + D})[2a](Di — Do) +(=D4/Dy)'), 
and 
¢=(r/DyD,)*(Di +Db)(D; —Dp). 


We recall that g and f’ must be such that dF/d6<0. 

It is evident that, despite these various restrictions, there still remains a 
wide range of permissible f functions. Figures 1 and 2 represent some simple 
exact solutions for the case D.=1, D,=9, f’=0, and several values of the para- 


meter a. 
We have here treated the case with both Dy and D, finite non-zero. Evidently, 
when only one of these quantities is required to be finite non-zero, the conditions 


on F are weaker and the problem is correspondingly simpler. 


VI. PracticaL UsE oF THE METHOD 

The typical mathematical problem arising in concentration-dependent 
diffusion is to compute ¢(8) when D(@) is known (frequently determined experi- 
mentally). Effective use of the method developed here depends then on the 
ability to match any given D(6) function with a member of the family of D(®) 
functions yielding exact solutions. 

This may be done, though rather prutally, by introducing (for the case of 
diffusion subject to (1.2)) an F function of the form 


F=2D; inverfo 0+ Ba,(1—O). ... seen ee (6.1) 
1 


D(8) funetion follows from (3.4). We may now mateh this 


The corresponding 
ers, Dy and a,, at (n+1) 


computed D(9) with its (n+1) adjustable paramet 
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points of the given D(6) curve. The problem reduces to solving (n ae simul- 
taneous equations for (n+1) unknowns. This is a possible, but not particularly 
attractive, manner of securing the solution. A similar method may be used 
for solving (1.1) subject to (1.3). 


o(N=s 


Fig. 1.—Parametric family of D(@) functions defined as follows : 
1>0> 3, D=9+ (1—6)[(2a+ 6 /2)0—3 / ]{2a+6+/r 12/B(4—40)}—3(2a+ 64/r)B(4—46) ; 
0<0<}, D=1—O[(2a+ 3v/7)0—(2a+ Vr) {20+ 34/—4/3B(40/3)}—3(2a-+ 3-V 7) B(40/3). 


Numbers on the curves denote values of a. 


It appears to the author that (both in the case of (1.1), (1.2) and of (1.1), 
(1.3)) itis preferable to keep the number of disposable constants in the F function 
down to perhaps two or three, and to provide the desired range of shapes of 
D(§) function by making use of the variety of the possible functional forms 
which F is free to assume. This requires the construction of a “ library ”’ in 
which are arrayed various F functions together with the corresponding D 


. ee ey re, 
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functions. It should then be possible, by inspection, to select an appropriate 


functional form of F, and the problem then reduces to matching the given D 
function by solving for the small number of disposable parameters occurring in 
the F function. It is beyond the scope of the present paper to provide such a 
library. Its construction and presentation remain a task for the future. 


Fig. 2.—Parametric family of ¢(6) functions defined as follows : 

1>0>3, 9=—6 inverfe (4—40) 4+a+4(a+3+/7)(6—2) ; 

0<0<?, p=2 inverfe 40/3+a+4(a+ 44/1)(8—32). 

Numbers on the curves denote values ofa. These curves represent 

exact solutions of (1.1), (1.3) corresponding to the D(§) functions 
of Figure I. 


VII. DISCUSSION 

In the present study we have made use of (4.1) and (5.1), Pi eet 

difficulties due to D(@) becoming zero or infinite at one or both ends a e 

9-interval ; but we have not provided a thorough investigation of such Fe a 
as the relationship (for diffusion subject to (1.2)) between the manner in whic 
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9(8) behaves as 0-0 and the behaviour of D(@) close to 6=0. It is hoped to 
return to questions of this nature in a later communication (cf. Philip 1957 for 


an early approach). 
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THE FUNCTION INVERFC 0 
By J. R. PHiir* 
[Manuscript received. November 11, 1959] 


Summary 


The function inverfe @ arises in certain diffusion problems when concentration is 
taken as an independent variable. It enters into a general method of exact solution 
of the concentration-dependent diffusion equation. An account is given of the properties. 
of this function, and of its derivatives and integrals. The function 


B(0)=(2/72) exp [—(inverfe 6)?] 
is intimately connected with the first integral of inverfe @ and with its derivatives. 
Tables of inverfe 6 and B(6) are given. 


I. INTRODUCTION 
The solution of one-dimensional diffusion problems is usually sought in the 


form 5 
concentration explicit function of distance and time. 


It has become increasingly evident, however, that there are occasions when it 
is simpler, and more illuminating, to seek the solution in the form 


distance explicit function of concentration and time. 


In particular, the latter approach has proved fruitful when applied to concentra- 
tion-dependent diffusion (Philip 1955) and when applied to problems where 
concentration-dependent diffusion is combined with a first-order (not; necessarily 
linear) phenomenon (Philip 1957). : 

In these connexions, it was found convenient (Philip 1955) to introduce the 
notation “ inverfe ” to denote the inverse of the function 


erfe onal exp (—C?)dG.  --- +225 ++ +e (1.1) 
ty ry 


Until now there has been no urgent need to examine in detail the properties of 
the inverfe function. However, inverfc 9 and its first derivative and first 
integral with respect to 6 enter intimately into the recently found general method 
of exact solution of the concentration-dependent diffusion equation (Philip 
1960). This account of the properties of inverfc 0, its derivatives, and its 


therefore forms an essential supplement to Philip (1960). The tabula- 


integrals 
oe here will frequently be needed 


tions of the functions inverfe 0 and B(8), given 
when the method of Philip (1960) is applied. 


* Division of Plant Industry, C.8.1L.R.0., Canberra. 
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The concept of inverfe as the inverse of erfc leads to definition of the function 
through (1.2) : 
@=erfe (inverfe 0). ......55-.-<2e5. (122) 


An alternative de novo definition of inverfc follows from equations (3.4) 
and (4.1) of Philip (1960). In this way the function may be introduced as the 
solution, F=inverfe 6, of the equation 


: _— 1 
om Spade EGS hsiasietne wares 1.3 
§ oe 29 ( ) 
subject to the conditions 
BPj=03 020218 a een eee (1.4) 


The following elementary results come directly from the known properties 
of erfc: 
inverfe 0=-+00; inverfe1—0; inverie2——o, .. (1.5) 


invvertc (2—0)=—inveric 02.2. .- se. 5) eee eee (1.6) 
We shall deal almost exclusively with the interval in 6, 0<0<1; it isa 


trivial matter to extend the results to the whole interval 0<6<2 by means of 
(1.6). Note that inverfe 6 is defined only within the latter interval. 


TABLE 1 
inverfe § COMPUTED FROM SERIES (2.4) 
inverfe 8 Computed from Series (2.4) inverteb 
6 Exact Value 
Three Terms -Four Terms Five Terms 
0-6 0-3706458 0-3707451 0-3707570 | 0-3708072 
0:7 0: 2724424 0- 2724557 0-2724566 =| 0- 2724627 
0:8 0-1791423 0-1791431 0-1791431 0-1791434 
0-9 0: 0888560 0: 0888560 0- 0888560 0- 0888560 


II. POWER SERIES FOR inverfe 6 
We introduce the power series connected with erfc: 


x 9d ee 


$n4(1 —erfe 2) =2 — sch aipi eal dale (2.1) 
Putting 
Om Fr 0) ever ayaa eee ++ (2.2) 
and writing « for inverfc 0, we have 
ce Con ae 


=o — 


3 +515 = art 419 ere ec : 7 hd Cale Teta (2.3) 
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Suppose now, that inverfe 6 (i.e. 7) may be expanded as a power series in #. 
Then we may formally establish this series by equating coefficients of powers of # 
on each side of (2.3). The result is 


. 7 127 
f = 193 5 7 
inverfc 0=3+ 33 +39? Cay ie 


4369 ., | 
22680 | fe ss 


ea (24) 


No simple general expression for the coefficients is apparent. It will be 
shown later in Section IV that (2.4) may be derived directly from the Taylor 
expansion of inverfe 6 about 96=1. The series of (2.3) is uniformly convergent. 
Presumably the series of (2.4) converges for | #|<}x}. It provides a useful 
means of calculating inverfc 0 in the neighbourhood of 6=1. Table 1 gives a 
comparison of the exact value of inverfc 0 with that computed from the first 
few terms of series (2.4). 


ELL AsyMPTOTIC FORMS OF inverfe 0, 6 SMALL 
For large values of « (=inverfc 0), we have the well-known asymptotic 
result : 


ee ae OU) Oe ee oe eee (3.1) 
TL 
(3.1) is equivalent to 
x2 —log 6 —} log x2’, 
which has the continued logarithmic form 
x2 —log 0—} log [x(—log 0 —} log [. . 
Accordingly we have the approximation for 9 small: 
inverfe 6={ —log 6 —4 log [x(—log eee aly Pork Le oe (3.2) 


As far as the author knows, no formal study has been made of the con- 
vergence of continued logarithms. The convergence of (3.2) is rapid for 6 small. 
See Table 2. In this table the symbol S,, denotes the nth member of the sequence 
formed by terminating the repeated logarithm at successive log 6’s. Thus, 


S,;=(—log 9)? ; 
S,=(—log 6 —4 log [7(—log 8)])* 5 
§,=(—log 0—4 log [x(—log 6 —} log {x(—log 9)})])*; 
and so on. 
It is evident that the limit to the accuracy of using (3.2) for § small is set 
by the limited accuracy of (3.1) rather than by me rave of convergence of the 
sequence S,. Note that S, proves a better approximation to inverfc @ than do 


the higher members of the sequence. 
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TABLE 2 
inverfe 8 COMPUTED FROM (3.2) 


Exact Value 


| | | we 

6 s S, Ss S, | inverfe 6 
| 
| 


NOY 3-7169 3-4540 3-4646 | 3-4642 3-4589 
Oss 3-0349 2-7437 2-7620 2-7608 | 2-7509 
Ome 2-1460 1-8081 1-8549 | 1-8480 . | 1-8214 


Pollack (1956) has established an inequality which leads to the following 
improved approximation 
2 exp (—2z?) 


“aa t/ (a2 +4]/n)} 6 eis le Oe wise wt ls 
This yields the better approximation 


6 


inverfe §6={ —log 6 —sinh-!V }n(—log 6—sinh—!4/[47(—log 6. . .}. 


(3.4) is more accurate than (3.2), and converges at about the same rate ; but 
it involves sinh-!, which is scarcely any simpler than inverfe. 


The following result follows from (3.2) or (3.3) 
ae inverfc 6 _ 
6—0 (—log Q)3 : 


IV. THE DERIVATIVES AND INTEGRALS OF inverfe 68 
Differentiating equation (1.2), we obtain 


a (inverfc 0)=—4x? exp [(inverfe-9)4]7 77... 00/2222. 2 ee (4.1) 


1 


ZT. ig convenient to introduce the function 


3) 

ray kd 
B(uj > 
TT 


exp. [== (ii veric:0)*)iw. | poeeieeoe eee (4.2) 


We then have 


cigs: Z. 
40 (inverfe Ves WEEK T eee Te, oe Sa (4.3) 


In addition, 


dB : 
"a0 =2 inverfc 0, 


and, in general (m0, n+0), 


d [(inverfe 6)™ 1 1 ‘ sh tahbeas | 
40 } Bn | =— Bazi l2n(inverte 0)"+14-m(inverfe 6)™—1]. et! (AS yepe 
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This result is useful for 2 te ! : 
generating the higher oh yA ; 
We have, for example, sg igher derivatives* of inverfe 0. 


al 


a aa 
ais (inverfc 6) = BR (2 inverfe 6), 
ahaa 1 
4 
f il : 
ae (inverfe Osa {48 (inverfe 0)?+-28 inverfce 0}, 
5 


der, 1 
ACE (inverfe 0) = — Be {384 (inverfe 0)4-++368 (inverfe 0)?+28}, 


do. i? 4 
FUL (inverfe 8) =F. {3840 (inverfe 0)5+5216 (inverfe 0)?+1016 inverfe 0}, 
oT (inverfe 6)=— 3 
qo? (inverte )=— Ri {46,080 (inverfe 0)®+81,792 (inverfe 6)* 
+27,840 (inverfe 0)?-+1016}. 
It. is evident that, for n odd, 
cag iL F 
aon (inverfe 0) = — Bn (a, —1 (inverfe 0)"—-!+a,_3 (inverfc 0)"—3+. . .+do), 
and, for n even, 
Ce ae ee 1 ; 
an (inverfe 0) = Bn {ay (inverfe 6)"—1 +a, _3 (inverfe 0)"—* 
+». . +a, inverfe 6}. 
In both cases the coefficients a,—1, etc. are all positive. 
Now 2/x!>B>0 throughout the interval 0 <§<2, whilst inverfe 0 is positive 
in 0<@<1, zero at §=1, and negative in 1<6<2. It therefore follows that, 
in the interval 0<0<1, 
aon (inverfe 6) is positive if m is even, negative if n is odd ; 
au 0=1, 
aoe, : ; : 
a0" (inverfe 0) is zero if m is even, negative if n is odd ; 
in the interval 1<0<2, 


n 
— (inverfe 0) is negative whether » is even or odd. 


* I am indebted to the referee for remarks which suggest the following, more elegant, treatment 


of the higher derivatives of inverfe 0. 


F=inverfc @ satisfies the equation 
d2F/d02=2F(AF/d0)®,  ..- ee eee e eee eee eee es (A) 


This may be established by differentiating (1.3) with respect to §. Now, if P,(f) denotes a 


a ae aMP/d0"=Ppy(F)(AF AO)", vee eeeeeee eee eee ees (B) 
it follows by differentiation and use of (A) that 
anti p/d0"+1 =P, 41(P)(aF/d8)y" 1, 
where ‘ 
Py pa(F)=(2nF+d/dF)Py(F). eevee ee eee eee (C) 


=], and P,(f)=1. Therefore, P,,(F) for all n>1 follows at once from 


Now (B) is true for 
g the higher derivatives of inverfe 0. 


(C), giving a simple means of determinin, 


B 
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We have the particular results : 


(—1yn (inverfe 0)= 0, 


do” 
Ch aa 
aon (inverfe 2) = — 0. 

The values of the first nine derivatives of inverfc 6 at 6=1 are 

—1int, 0, —2(4n4)3, 0, —28(474)5, 0, —1016(47:4)’, 0, —69,904(3704)°. 
It follows from these results that equation (2.4) may be established by applying 
Taylor’s theorem to the right-hand side of the identity 
inverfe §—inverfe [1 —(1—9)]. 
We also note that it follows from (4.4) that 


d"B dz —1 


pale 4) ee | ee 4.6 
aon 29a (inverfe 6) (4.6) 


It is readily established by integration by parts, or by use of (4.1) in (1.3), 
or by integrating (4.4), that 


i) 
| inverte 0 ao=- exp (—inverfc 6)?=3B. sing ASRS 
0 
A further integration yields 
if i inverfc 6 d§d0 =—_——. Taw = erfe (4/2 inverfe 6). 


This and the higher integrals of inverfc 6 do not appear to be of significance in 
the present developments. 


V. THE Function B(8) 

We have seen that B(@) is simply related to the first derivative, and to the 
first integral, of inverfe @. For this reason it proves of primary importance 
in the development of the general method of exact solution of the concentration- 
dependent diffusion equation (Philip 1960). 


We note from (1.6) and (4.2) that 
B(2—0)=B(6). W518 ote hereto (5.1) 


We have already remarked on the simple relation between derivatives of 
inverfe 0 and those of B. It follows that the Taylor expansion of B(8) about 
§=1 yields as the power series for B 


2 4 7 127 4369 
B(§)= yon Ae Cie er igh 
(9) if o 6 90° — 520° — 773400" .4 | ap (Osa) 
where # is Genin defined by (2.2). Presumably the series of (5.2) converges 


for | #|<474. It enables B to be calculated readily in the neighbourhood of 
O=1, 
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The behaviour of B near 0=0 is of interest. Approximation (3.1) applied 
in (4.2) yields 
BO) 220 log 6) 4pians) eee iow artes (5.3) 
and it may be shown that 
B(0) 
lim —————_.\ = 1. 
iouae ao (—log 0)? : 


It follows that, as 6+0, B—-0 more rapidly* than does §1—*, where ¢ is any non- 
zero positive quantity, and more slowly than does 0. 


VI. TABLES OF inverfe 0 AND B(@) 
The only existing table of inverfc § known to the author is in Fowle (1921, 
p. 60). The column of the table headed v/e gives @ in the present notation, 
and that headed 2q gives values of 2 inverfe 0 to four or five significant figures. 
In connexion with Philip (1960) it is helpful to have a table of inverfe @ readily 
available. In the course of constructing the table of B(0), it was a simple matter 


TABLE 3 
THE FUNCTIONS inverfe 80 anp B(@) 


| \ 
) inverfe 6 B(0) 6 | inverfe 8 B(8) 
0 oo 0) 0-40 0595) L167 IS) 07915 185159 
LOSS 3:458 911 0-000 007 177 8 0-45 0-534 159 1 | 0-848 280 6 
10-6 3-122 587 0-000 065 742 0-50 0-476 936 3 | 0-898 807 9 
MO 2-750 936 0-000 583 27 0-55 0:422 680 2] 0-943 766 6 
Ose 2°326 754 0-005 026 6 0:60 0-370 807 2) 0-983 423 2 
Os 1-821 386 0-040 898 3 0-65 0-320 858 31] 1-017 992 1 
0:05 1-385 904 01655 219 95 0-70 0-272 462 71] 1:047 646 6 
0-10 1-163 087 02291 S71 6 0-75 0: 2255312 ee lcO72ne oc Giant 
0-15 1-017 902 0:400 379 2 0-80 0-179 143 4] 1-092 741 7 
0-20 0-906 193 8 | 0-496 384 2 0-85 OSS Se 2O eon be LOSms 1Omms 
0-25 0-813 419 8 | 0-582 241 7 0-90 0:088 856 07] 1-119 505 3 
0-30 | 0-732 869 1 | 0-659 472 6 0-95 | 0-044 341 3] 1-126 162 8 
0-35 | 0-660 854 4| 0-729 098 6 1-00 0 1-128 379 2 


to develop a new and more accurate table of inverfc 0. Details are given below, 
and the resulting tabulation is presented in Table 3. No graph of inverfe @ is 
given, since the shape of erfc a is well known. , 
No table of B(0) is known to the author. The tabulation of this function, 
which is also given in Table 3, was constructed with the ald of National Bureau 
of Standards (1954) tables by methods described below. Figure 1 gives the plot 
of B(0). 
Table of inverfc ®. The table of inverfe § was constructed from the National 
Bureau of Standards (1954) tables of erf « by a process of inverse interpolation. 


Eeuppoe a FEO ANB) =) Thon wo tay EAU ga ERS ISNT 
09 8-0 


P 
than does Q, provided that , oe —(0)=0. 
)-—+00 
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“ Linear ” interpolation of the type suggested in the introduction to those tables 
proved sufficiently accurate to ensure that errors in the final place given would 
not exceed unity. Twenty comparisons with Fowle’s table were possible ; in 
every case the final place of Fowle’s tabulation of 2q (i.e. 2 inverfe 6) was 
confirmed. 


o-8 


0-4 


0-2 


H——_> 


Fig. 1.—The function B(6). 


Table of B(®). Once inverfe 0 was computed, it was a simple matter to 
calculate B(0) by linear interpolation (again of the type suggested in the intro- 
duction) in the tables of the derivative of erf «in the National Bureau of Standards 
tables. It was established that this process would not yield errors greater than 
unity in the final places shown in the table. 


Most of each table was checked by differencing. 
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THE DISTRIBUTION OF IONS FORMED BY ATTACHMENT OF 
ELECTRONS MOVING IN A STEADY STATE OF MOTION 
THROUGH A GAS 


By C. A. Hugst* and i. G. H. Huxipy+ 
[Manuscript received December 3, 1959] 


Summary 
The distribution of ions formed by attachment of electrons diffusing through a gas 
is solved exactly, and the results compared with an approximate calculation given 
earlier by Huxley. The corrections to the approximate results are inside the present 
experimental error, and so confirm the satisfactory agreement with experiment already 
obtained. 


I. INTRODUCTION 

The distribution of electrons and ions in a stream from a pole source and 
drifting and diffusing in a gas when the processes of ionization by collision and 
attachment to molecules are active, has been discussed in another paper (Huxley 
1959), where an exact solution is given for the distribution of the electrons 
and an approximate solution for that of the ions. It was apparent that under 
conditions where the mean energy of agitation of an electron was as great as 
20 times that of a gas molecule, or greater, the approximation to the exact 
solution was close and that the theory could provide the basis of an accurate 
method for measuring the coefficients of attachment of electrons as a function of 
mean energy of agitation in a given gas. 

This expectation has been justified by the application of the method to 
oxygen in which self-consistent and accurate measurements of attachment 
coefficients were made (Huxley, Crompton, and Bagot 1959). 

It is important, however, to extend the measurements to the range of small 
energies of agitation of the electrons from approximate thermal equilibrium 
to about 20 times this value, and for this purpose an exact solution of the problem 
is required. This solution is given in what follows. 


II. DISTRIBUTION OF ELECTRONS 
Consider a pole source of electrons placed at the origin of coordinates in 
a gas and emitting electrons at constant rate of S per second. Let W be the 
drift speed of electrons through the gas in a constant and uniform electric field # 
and let D be the coefficient of diffusion of the electrons. Suppose that the 
of ionization by collision and electron attachment are both operative 


@88es 
aa If the positive direction of the axis Oz 


with coefficients «, and «, respectively. 
* Department of Mathematical Physics, University of Adelaide. 
+ Department of Physics, University of Adelaide. 
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be chosen to be that of W, the differential equation satisfied by concentration n 
of electrons in the stream is (Huxley 1959) 


V2m=2AOn/O2z+2r00N, «ee eee ee eee eee (1) 
where 
2A=W/D; a=a,—«;. 


The solution of equation (1) appropriate to the pole source at the origin is 
(Huxley 1959) 

n(@,y,2)=(S/47-Dr) exp (A@—ur), = « - + 0 ee eee (2) 
where 

uta=,24-205 fa y*te*. 


The distribution n’ due to a pole source of strength S exp (2h) placed at the 
point (0,0,2h) is 


n' (x,y,z) ={S exp (2Ah)/4xDr’} exp {A(2—2h) —ur’} 
==(S/4Dr') exp (A ur), cece en eee eee (3) 
where 
vr ={a? +y? +(¢—2h)*3, 


It follows from equations (2) and (3) that the solution of equation (1) that 
represents the distribution n(z,y,z) in a stream of electrons emitted at the rate 8 
from a pole source at the origin and drifting and diffusing to a metal electrode 
coinciding with the plane z=h, over which n is now everywhere equal to zero, 
is (Huxley 1959) 


n=(S/4rD) exp (Az): {exp (—ur)/r—exp (—ur’)/r’},  2<h, 


since r=?’ over the plane z=h. 


III. DIsTRIBUTION OF IONS 
Consider the distribution of negative ions formed by attachment of electrons 
to molecules in the stream of electrons from the isolated pole source represented 


by equation (2). Let N(a,y,z) be the concentration of the negative ions. The 
differential equation satisfied by N is, 


V*NS=2),0N {02 —(o Wald. ane ee (5) 
where 


2A, = W,/D,, 


and W, and D, are the drift speed and the coefficient of diffusion respectively 
of the negative ions. The ratio A/A, is approximately (or, with a Maxwellian 
distribution of agitational speeds, exactly) equal to the ratio of the energies of 
agitation of ions and electrons respectively. 


It is required to find the particular solution of equation (5) given that n is 
a solution of equation (1). 
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Put 
z 
P=C exp (—#e) { OXD(h2.) 202 « Batemans eee (6) 


in which C and &k are constants. 
It follows that 
oP /dz= —kP +Cn, 


and | Ne eae se a cers (7) 
@P/de2—k2P —Okn+Can|az. | 


Also, consider 


— oO 


I z 3 O2n 2 
=C exp (—kz) exp (kz Jade Mb dude eo ete (8) 


and let exp(kz’)én/dz and exp(ke’)n both approach zere as 2’ tends to —oo. 
After two successive integrations by parts and comparison with equation (7) 
it can be seen that 


l=¢*P/o2*, 
and also 


V?2P=C exp (—He) | exp tke’) Vilnideens: wt, team (9) 
To solve equation (5) substitute the value N =P in equation (5) and use equations 
(7), (8), (9), and (1), to find 


of” exp (kz’){220n/dz’ +2ran +2, kn}de’ +(%,W/D; —2,C)n exp (kz) =0. 


On differentiating this equation with respect to 2 and collecting terms, it reduces 


to 
(2ra0 +a,Wk/D,)n +{a,W/D,—2(A, —))C}an/dz=0, 


whence 
ie iW 3 AG, LY 
2004 —A)D, (pA) D,’ agate. ateveteieters (10) 
k= —)Aa/(Ay—A) 

Thus 


Aa D Aa # ( Aa ) dz’ 11 
aa a 2 exp | — 2’ \ndz’. .. (11) 
a, Dee ee Nhs a eee 


When the solution of equation (1), n—S exp (Az—ur)/4Dr, is adopted, then 


Sar Ag ‘bg ex es sel rae ee 
ee eee Les A p —wU 
Se ae ren : i\ 
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The distribution of electrons from a source at the origin and its image at (0,0,2h) 
is given in equation (4) and the corresponding distribution of ions is 


Sa,A ha z f o : 
weer eee ieee exp i — ——=4%e 
De es (esa P | ae 
< {exp i) exp Ce Nae, z<h 
5 


where r’ =(o2-+2'2)3, r”={p?-+(2—2h)?}#, and p? =a? +y?. 

The current density J at a point (p,h) on the receiving electrode is 
J=(eNW,),-, and when A,>A this expression becomes essentially the same as 
that derived by Huxley (1959, equation (25)). 


An alternative solution of equation (5), namely, 


Sia Ae Ne i x é 
=-, ea oe a p Remreeny Fe 
USTED eee eerie "ae ( 3) 
- =P ee — N de’ 


is inadmissible since in it N does not vanish at z= — oo. 

Equation (13) does not make N=0 at z=h, and so does not represent 
the exact solution. If a suitable solution of the homogeneous equation corres- 
ponding to (5) is added then this boundary condition can be satisfied. The 
resulting expression is rather complicated and can be evaluated only approxi- 
mately. The physical basis of the approximation may be represented as follows : 
suppose a uniform stream to approach a plane electrode at right angles to the 
stream, the concentration at a large distance from the electrode being N,. 

At a distance s from the electrode, the concentation is 


N=WN,{1 —exp (—2),s)} 


and the flux of ions to an area dS of the electrode is 


(7.53) dS =2),D,N,d8=W,N,d8, 
s=0 


which is the flux across an element dS of a geometrical plane normal to the 
undisturbed stream. 

In an ionic stream 21, =W,/D,~40H, where £ is the electron field strength 
in volt cm~' and W, and D, are expressed respectively in em sec! and em? see, 

If H=1 V cm-1, it follows that N~2N, at a distance of 0-5 mm from the 
electrode. 

It follows that the distribution in the stream of ions of equation (13) in 
the vicinity of the electrode at z=h is given closely by N {1 —exp —2A,(2h —z)} 
and that the flux to the surface element dS of the electrode at a point (p,h) is 
NW,dS, where WN is given by equation (13) with z=h. This, as explained 
above, is the procedure that was adopted in practice. 
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The exact mathematical treatment is most conveniently expressed in terms 
of Fourier Bessel transforms. The approximate solution (13) can be written 
in the alternative form 


N _ «,WS ae ce) hdl J (ke) \ exp (¢—2h)4/(y2+k?) 
Be ils VA RTOS Ate (uk) 
exp {—2v/(w?+k*)} ) (13’) 


aA+(A—A,){A—4/ (we? +k)} ‘a Sanco. dO, Cone OLG.0 Boat © 


The expression (13’) is defined for z>0. For z<0, the second term inside 
the bracket must be replaced by exp {24/(u?-+k?)}/[aA+(A—A,){A + +/(u? +h?)} J. 
The two parts of the solution so obtained do not join up smoothly at z=0. 
But suitable choice of the complementary function for <0 will ensure that the 


solution is continuous and satisfies the correct boundary condition at z=—o. 
So strictly speaking (13’) is not equivalent to (13) but rather to 
= Sa,A BALE \ Nas, % , exp (—yur") 
as a cp | —— }- d a 1X7 ; 
Same), (4) | We ira ( a r 


eae le eee ly ee 
is dz exp i oa w \ 


It can be verified that the following choice of complementary function is suitable : 


“WS (# Jo(ke) 
= & k Ls zs, 
SDD, Vet 


xexp [(2—h) {Ay +-V/ (Ar +h*)} +h{a— Vw? +h)} ] 


re 


1 i 
| aOR VET Ree | 


RO ray ae: (14) 
The ion current at the anode is now 
—D, LW +¥,) | 
102 een 
From (13), 
oN AaD, 
a Vi ; 
a 02 \e=h Ay—A z=h 


The calculation of @N,/d2 |--» can be carried a approximately if > or 
. Ap. The integrand contains a factor A; ++/(Ai+4?) on diflerenGaning: ek 
respect to z, and then putting 2=h. If A,>w, we Bele Ay gf Vv shoae 
= + V/{(At—p2) +(u2+4*)}, and expand in powers of (y2+4?)/(Ai—p?). To 
the first order we have 

2 Pe oe VO ee 
mt VIM+HWt VIM HOF oe ia) VOR — we Va +) + V2} 
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Neglecting the last term, this gives the approximate expression 


ON 2 / dD, Cres | 
—D,=— ==) Ai —2)} —,——— ——e ”N 
D, Oz — {4a + v/( 1—U )} rie /02—p?) a2 we 
~2n,D,N | +71(sae ew ee ) .. (15) 
lz=h dz lz=h lz=h 
agreeing with the physical argument given earlier. 
If A;~p, we write 
1 —p2 
Art (At +h?) =A + (U2? +h) ta oa) (uz +I) 
(At —p2)? 1 


TIVE TR) 6/4) + VTP 


This is an expansion in powers of (At —u2)/(u2 +k), and, to the first order, we 
can neglect the last term. Then, with the notation introduced above, 


ON |", | sei aT 
=D,| a )_ =npw DS a NaN 


oe} h 
+4D,03—u2)( | dz’N+— | az'N-). 
h 


The expressions (15), (15’) can be calculated numerically without too muct 
difficulty. 
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Il. THE INITIAL DESIGN OF COMPOUND TRIPLET SYSTEMS 
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Summary 


Part II. The general method developed for the type 111 triplet is extended to 
the initial design of triplet systems having compound components. The possibility of 
replacing a single lens by an equivalent doublet or triplet allows the single lens to be 
regarded as not being restricted to values of N and V associated with known glass 
types. It follows that the glass parameters «, 8, y, and — may be treated as continuously 
variable over limited ranges. The general effects of the variation of these parameters on 
the powers and separations of the basic initial solutions are shown in a set of diagrams. 
As an example, the initial design of a triplet with compound members is given in detail. 


I. INTRODUCTION 

In Part I of this series (Cruickshank 1958) the basic equations and procedures 
have been given for the initial design of the parent triplet objective, type aie Te 
for a given selection of glasses and prescribed values of the power, the paraxial 
chromatic aberrations, and the third-order monochromatic aberrations. The 
purpose of the present paper is to extend the general method to the initial design 
of other members of the triplet family (Fig. 1 of Part I) having compound 
components. The first step necessary is to consider the equivalent glass properties 
of such components. 


II. EQUIVALENT DOUBLETS AND TRIPLETS 
Consider the thin lens (9,N,V) of power ¢ made from glass of refractive 
index N and V-number V. Suppose this lens is to be replaced by n thin lenses 
(ky9,N1,V 3), (kop,No,V2), - + +5 (k,9,N,,V,) in contact, the replacement being 
made in such a way that the compound group has the same power, the same 
Petzval sum, and the same paraxial chromatic aberrations. This requires 


that 


» k; eer eR RMN NC uci (1) 
j=l 
k;@;=@, Sie prada ie! <a" 0: regen ener a: Sue sige (2) 
g=1 
Sh Weds ee eae ec enina c (3) 


where w=1/N, ,=1/N,, b=1/V, p=1/V; Two particular cases arise for 
detailed consideration. 
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The first of these is the case in which the replacing group is a doublet, i.e. 
the thin lens (9,N,V) is replaced by the thin lenses (k,p,N1,V;), (kop,N o,V 5) 
in contact. The parameters k,,k, must be chosen so that 


kg-bhig a1, Pak ee eee eee (4) 
and either 
Cyl, EO ghg = Oy 9 tne dee ee eee ie (5) 
or 
Giga, HO en eae nee eee (6) 


Many possibilities now arise. Suppose, for example, that the single positive 
thin lens (9,NV,V) is replaced by a doublet consisting of a thin positive lens 
(k,o,N,,V,) and a thin negative lens (k.9,N.,V.) in contact. Then k, will be 
negative, | k./k, | <1, and both V, and NV, must be greater or less than N according 
as N,2N,. Similarly, both V, and V, must be greater or less than V according 
as Vso Vi. 

Looking at the matter from the opposite point of view, suppose a thin 
doublet of positive power is constructed from the positive thin lens (k,9,N,,V,) 
and the negative thin lens (k,9,N,,V.) in contact. This doublet is equivalent 
as regards power, Petzval sum, and paraxial chromatic aberrations to a single 
thin lens (9,NV,V) where N and V are given through equations (5), (6). These 
values V,V will not correspond in general to those of any known glass and may be 
varied widely by intelligent choice of the two glasses and the value of k,/k,. 
In particular, by selecting a pair of glasses of the same index for some mean 
wavelength and varying k,/k,, the effective value of V may be varied over a wide 
range without change of index. Such an iso-N doublet can be very useful. 
For example, if we combine a positive thin lens of the DBC glass 658508 with a 
negative thin lens of EDF 654335 then the equivalent thin lens has an effective 
N-value of 1-66 approximately and an effective V-value anywhere between 50 
and infinity depending on the value of k,/k,. Similarly, by choosing an iso-V 
pair of glasses the effective value of N may be varied substantially without 
change of V. More generally, the choice of glass pairs in which both NV and V 
are dissimilar allows the effective values of N and V to be changed simultaneously 
over a considerable range. 


Consider next the special case in which the replacing group is a triplet. 
Thus suppose the thin lens (9,N,V) is replaced by a compound triplet consisting 
of the three thin lenses (k,9,N,,V,), (kop,N2,V2), (kg9,Ng,V3) in contact. The 
replacement may now be made in such a way that the triplet has the same power, 


the same Petzval sum, and the same chromatic aberrations as the singlet it 
replaces. Hence the equations (1)-(3) become 


kythk.+k, . =P OME ene ie me (7) 
Wy hy -F alts + Gighg== G08 oer ara ie eee al (8) 
Pky + Vals Uskg =U," Geeeeme eee es (9) 


whence 


hy =AAs kyeSAiAs Sk Ac rrr (10) 
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where 
of al i i 1 et 
ee! Or Oa @s |, AL=| 0, @ | @s |, 
vb de de Ids be | 
a 1 il 1 il! dlp 
ies WaiaeGe Gu. |, S| Foy A Oye 
it, te ¢ id, te Gs 


It is thus formally possible to design a triplet component which is equivalent 
as regards power, Petzval sum, and chromatic aberrations to a thin lens having, 
within a certain range, any value of NV, and, independently thereof, any value 
of V. Therefore in the development of a photographic objective we may ignore 
the restrictions imposed by the existence of only a limited number of glass types 
if we are prepared to use cemented doublets and triplets to replace components 
which are single thin lenses in the initial stages of the design. For example, 
if for some reason one requires for the negative component of a triplet a glass for 
which N=2-0 and V=15, it can be assumed that this fictitious glass exists and 
in the next stage of design these effective values can be achieved by replacing 
the fictitious component by an actual cemented triplet. : 

In the design of triplet objectives of type 111 we are limited to available 
glass types for the three components and hence the values of the glass parameters, 


a=V,/V,, Ses 
P=N,/N,, y=N IN, » 


are not continuously variable quantities. In the light of the considerations of 
this section, however, those members of the triplet family of objectives which 
have compound components may be considered at first as simple type 111 triplets 
in which one or more components have fictitious glass constants. This simplifies 
the problem of designing such systems as the theory of the initial design of 
type 111 triplets given in Part I is immediately applicable since each equivalent 
simple triplet has the same power, the same corrector power y, and the same 
Petzval sum and chromatic aberrations as the compound system to which it is 
equivalent. The appearance in this way of glass constants which do not 
correspond necessarily to known types removes substantially the restrictions 
on the glass parameters and it is neither artificial nor without physical meaning 
to regard a, B, &, and y as continuously variable over limited ranges at least. 


TII. Tue EFFECTS OF THE VARIATION OF THE GLASS PARAMETERS 

The first stage in the development ot the design of a triplet objective is 
the determination of the powers 9,, 2», ?, and the separations t,, t, of the three 
thin lenses a, b, ¢ constituting the initial arrangement of the system. From the 
point of view of this first stage the variety of constructions in the triplet family is 
largely the consequence of the variation of the glass parameters CH. é, and Y. 
If an understanding of the range of properties of the whole family of objectives 
is to be reached, or, indeed, an understanding of the detailed properties of any 
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one member, it is necessary to know the effect on the initial solution of the 
variation of a, 8, &, and y. It is not assumed that this is the whole story as 
regards glass selection, because each of these parameters is a ratio and therefore 
any given values of the parameters might be achieved in a variety of ways. 
However, it represents a systematic approach to what is a complicated problem, 
and, once this general approach is outlined, it may be supplemented by a closer 
study of the effects of achieving any of these parameters in different ways in 
connexion with different members of the triplet family. 


When the changes in the powers and separations due to given changes in the 
glass parameters are known there remains still the problem of determining whether 
the new basic solution is more favourable or less favourable than its predecessor. 
This is effectively the second stage of the design. Broadly speaking, an initial 
solution must be scrutinized in three different ways. The inherent characteristics 
of the system in respect of spherical aberration of all orders must be examined, 
then the asymmetry errors, i.e. the comas of all types and orders thereof, and 
finally those characteristics which control the field correction, i.e. the astigmatism 
of different orders and what is generally of over-riding importance in photographic 
objectives the oblique spherical aberrations. 

In this laboratory extensive use is made of the aberration coefficients 
introduced by Buchdahl (1954). The general method developed for the use of 
these coefficients is described in another paper by Cruickshank and Hills (1960). 
Anyone who has used these coefficients will readily concede that they provide 
the most powerful means available at present in geometrical optics for predicting 
the image quality of a system and for analysing the way in which the particular 
construction achieves its performance. The laboratory’s computing room 
furnishes a standard service of the computation of all the aberration coefficients 
of third, fifth, and seventh orders for any system. Formerly these coefficients 
were computed on electric desk calculating machines, but recently Ford (1959) 
has arranged the computation on a digital computer. In these computations 
the total value of each aberration coefficient of the complete system is obtained 
as the sum of a set of partial coefficients, one for each surface of the system. 
If, for example, a variation of the basic parameters leads to a general reduction 
of the partial coefficients of the spherical aberration at the several surfaces of the 
system, it is permissible to infer that this represents a step towards the develop- 
ment of reduced zonal spherical aberration and hence the possibility of increased 
aperture for the system if needed. One has only to compare the partial aberration 
coefficients of a few first class objectives designed for apertures of, say f/3-5, 
f/2-5, f/2, and f/1-5 respectively, to see that this is so. No successful high 
aperture systems are developed from initial arrangements with characteristically 
large partial coefficients of spherical aberration. Now, broadly speaking, 
reduced partial coefficients of spherical aberration go hand in hand with reduced 
powers of the components, which in turn lead to shallower curves throughout 
the system. Of course, reduced powers cannot be sought without consideration 
of the separations of the components which may become so small that the system 
cannot be built or so large that there is an unacceptable degree of vignetting. 
The criteria by which the inherent properties of a system in respect of asymmetry 


Be 
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errors and field errors are to be judged are Similarly based on an examination 
of the appropriate aberration coefficients. The calculation of the aberration 
coefficients provides an analysis of how these aberrations arise within any 
particular system, but just how these are linked with simpler quantities in the 
initial solution requires deeper investigation. Further work is in progress in this 
connexion, in the course of which it is planned to make a broad study of the 
characteristic distribution of these aberration coefficients within the various 
types of construction which constitute the triplet family. 


(a) Objectives for which €=y=1 

Consider first those objectives in which the effective glass constants of the 
front and back components are the same. It is common practice to construct 
type 111 objectives from two glasses, using the same crown glass for both positive 
lenses. This makes £=1=—y and leaves only « and £ as the effective glass 
parameters. Extending this idea to triplets with compound components it is 
possible to make systems like types 212, 222, 313, 323, etc. on this pattern 
provided the effective N and V values for the front and back components are 
the same. 

Choosing some typical values for the aberration residnals which enter into 
the fundamental equations 

X= —0-25, L=-+0-20, 

P=-+0-50, T=-+0-05, 
the variation of the powers and separations of the thin lenses in the initial 
arrangement have been calculated for a range of values of « and 6 that are quite 
practical. Figures 1 (a)-1(f) present a summary of the results. 

It is immediately clear that increasing the value of f£ leads to substantial 
reductions in the powers of the first two lenses and to a lesser degree reduces the 
power of the back lens also. Intelligent choice of 6 can lead, therefore, to 
shallower curves throughout the system, with consequent reductions in the 
values of the partial aberration coefficients at the individual surfaces of the 
system. As 8 increases from lower values the front airspace at first widens but 
subsequently diminishes, while the rear airspace decreases steadily throughout 
the range. For any value of @ the | 9,| versus « curve exhibits a shallow 
minimum somewhere towards the higher end of the « range. In addition, 9, 
is decreasing steadily as « increases, while the accompanying increase in 9, 
is not nearly so large. This may suggest, in general, the selection of an a-value 
near to that at which | 9,| attains its minimum. However, the very consider- 
able increases in t, and f, as the value of « increases must be weighed acest this. 

There is another important matter to be considered in connexion with the 
choice of «, and that is the correction of the spherical aberration of the system. 
It was shown in Part I of this study that, once the glass types and values of 
desired residuals have been chosen for a triplet, the shapes of the three components 
can be adjusted to secure control of the primary astigmatism, distortion, aut 
coma, while the primary spherical aberration, oj, 18 controlled by the parameter 
the power of the corrector of the system. The whole placing of the o, versus x 


parabola is often controlled by the choice of « as is seen in Figure 2. 
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It cannot be assumed, however, that everything is necessarily ne 
provided the co, versus y parabola comes near to or intersects poe As AS « 
decreases and the corresponding parabola drops towards the — the joe 
magnitudes of the partial aberration coefficients of the surfaces increase. Varia- 
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Fig. 1.—The effects of the variation of 6 on the powers and separations of the thin com- 

ponents of unit power triplets having different values of «~ are shown. In all triplets 

&=1, y=1, and X=—0:25. The residuals of Petzval curvature and the chromatic 
aberrations are in all cases P=0-50, L=0-20, T~0-05. 
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tion of x for a fixed «-value leaves the partial coefficients of spherical aberration 
of the negative lens almost unchanged, while the partial coefficients of the positive 
lenses change considerably in accordance with the parabolic law. On the other 
hand, the partial coefficients for the negative lens, and hence also those for the 
positive lenses, increase very substantially as « is decreased. This is due 
principally to the effect of « upon the front airspace. A decrease in « lowers the 
front airspace, raises y,, and hence the spherical aberration coefficients at the 


PRIMARY SPHERICAL ABERRATION 


POWER OF CORRECTOR, X 


Fig. 2.—This figure shows the effect of the variation of « on 
the placing of the o, v. y parabola in a type 212 triplet. 


surface of the negative lens. If, then, it is desired to secure as high an aperture 
as is consistent with a certain type of triplet construction, it is important to 
consider the location of the co, versus 7 parabola during selection of glass go: the 
system, for the lower order partial aberration coefficients largely determine the 
magnitude of the higher order aberration coefficients of the system. A thorough 
study of Figures 1 and 2 will enable the reader to grasp pone of the essential 
significance of glass selection in the design of this class of triplets. 


(b) Objectives for which Exy. The General Case 


Apart from the construction of a series of three-dimensional models, it is 


difficult to present the effects of the simultaneous variation of the four glass 
parameters on the powers and separations in the general case of three-glass 


; i to set B=1=y and 
has been done, therefore, is firstly 
oe fon of &. The results are presented for selected 


It is seen that an increase in the value of @ 
increases the powers of the lenses in the corrector system and eae aren 
5 changes being almost linear over this 
er of the back component, the ¢ 
sa A each case. Over against this are fairly rapid decreases in both airspaces. 
Secale setting B=1—&, the variation of y has been investigated over the 
? 
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investigate the effect of the varia 
values of « in Figures 3 (a)-3 (f). 
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range 0:8-1:2. The results are presented in Figures 4 (a)-4(f), from which 
it is seen that an increase of y increases all powers, decreases the front airspace, 
and increases the back airspace. At smaller values of « the effect on the front 
airspace is very small. It is noticeable in all these variations that the effects 
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thin components of unit power triplets having different values of «. In each triplet B=1, y=1 
and y=—0:25. The residuals of Petzval curvature and chromatic aberration are P=0-50, 


4 
(e) (f) | 
Fig. 3.—These figures show the effects of the variation of — on the powers and separations of the 
L=0-20, T=0-05. 
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on the powers of the first two lenses are very much greater than on the power of 
the back component. This component behaves as a fairly stable feature of the 
construction, while the front portion of the system may undergo wide changes. 
This strengthens the view adopted by the writer that the first two lenses fulfil 
the function of a corrector system. 
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Fig. 4.These figures show the effect of the variation of y on the powers and 
.. a of the thin components of unit power triplets having different values 
re each triplet B=1, 6=1 and y=—0:25. The residuals selected are 
. n > , 
“fa P=0-50, L=0-20, and T7=0-05. 
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IV. AN EXAMPLE OF INITIAL TRIPLET DESIGN 

It will be useful now to discuss in detail an actual example of an initial 
design to illustrate the general method. Suppose we consider what could be 
achieved in a type 212 lens in which the front and rear components are iso-index 
doublets. We choose from the glass catalogue the approximately iso-index 
pair of glasses 1-65695, 50-80 and 1-6538, 33-55 for this purpose. The equivalent 
refractive index of a doublet from this pair will be approximately 1-66 but the 
equivalent V-number will depend on the relative values of k, and k,. Equations 
(4), (6) show, for instance, that for the range 0<| k,/k, | <<0-6 the corresponding 
equivalent V-number range is 50-8<V<222, so that we may give V almost 
any value we please. 

Figures 1 (a)-1 (ce) show that the powers of all components are reduced by 
making £ as high as possible. In addition, the powers of the first two components 
are reduced by using a high value of «. Having fixed the value of V, at approxi- 
mately 1:66 by the choice of the iso-index pair, there remains only one way of 
increasing 6 and that is the choice of a low index glass for the central negative 
Jens. Accordingly we choose the BSC glass 1-5095, 64-4 for this. It may seem 
rather unusual to choose such a crown glass for the main dispersive member 
of the triplet, but at least it does provide the increased value of & which seems 
desirable. We do not know yet what value of « will be required to make the 
correction of spherical aberration possible. Initially, a high value of « is desirable 
to give further reduction in the powers of the first two components. Quite 
arbitrarily we choose «=2 and see what happens. Now «=V,/V, whence 
V,=2 X64-4=128-8. Equations (4), (6) therefore give k,/k,=—0-54083 and 
k,=2-17785, k,=—1-17785. With these k values equation (5) gives 
N,=1:66068 and hence B=N,/N,=1-1002. 

Since we are using the same iso-index pair for the back doublet, the equivalent 
refractive index of this doublet will be also approximately 1-66 and hence 
y=N,/N, will be very close to unity. It only remains, then, to settle the value 
of € Figures 1 (e) and 1(f), which are drawn for E=1, y=1, show that the 
choice of a=2 and B=1-100 results in high values for the airspaces if E=1. 
Figures 3 (e) and 3 (f) show, however, that if € is increased the airspaces will be 
reduced substantially. Figures 3 (a)-3 (c) show also that this reduction in the 
airspaces will be obtained at the expense of some increase in the powers of the 
first and second components. A compromise is therefore called for and so, 
tentatively, we put €=1:5. Since €=V,/V,, this gives V.—128-8/1-5—85 -867 
and thus k,=1-79439, k,= —0-79439, and k,/k,——0-44270. It then follows 
that the equivalent index of the back doublet is WN -=1-:65946 and thus 
y= NN = 150007. 

The next matter is to select the values of P, Z, and T. Putting P=0-6 
makes the Petzval sum P/N,=0-6/1-66068=0-3613. This is fairly high but can 
be tolerated if the field angle required is not too great. If we make L slightly 
positive the longitudinal chromatic aberration will be corrected at an outer zone 
of the aperture rather than paraxially. We therefore put L=0-50, and similarly 
T=0-045. Now, using equations (21)-(32) of Part I, initial solutions are 
made for three values of the corrector power, say y= —0-25, —0-375, and —0-50. 


Similarly, the con 
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The resulting equivalent thin solutions are summarized in the first five lines of 
Table 1. Using the values k,/k,=—0-54083 and —0-44270 for the front and 
back doublets respectively, the powers of the members of these doublets are 
calculated in each case and are shown in the next four lines of Table 1. This 
completes the calculation of the powers and separations of the thin components 
for these three systems. 


TABLE | 
DETAILS OF DESIGN EXAMPLE DURING THIN LENS STAGE 


x —0+250 | —0'375 |  —0-500 
4 1-2723 1-2016 1-1322 
2, —2-1306 _2-1710 —2+2156 
®, 1-6705 1+ 7857 1: 9040 
a 0-2244 02279 02326 
a 0-1351 0-1169 +1023 
Pa 2-7708 2-6168 2-4657 
Das —1-4985 —1-4152 —1+3335 
Qa —1-3270 —1-4184 —1-5124 
Peo 2-9975 3-2041 3-4164 
ome 01853 0+ 2276 02718 
S.. | 0-4992 0:4213 0-3401 
oe —0-2744 —0+ 2027 —0+1232 
S —1-9954 —1-9220 —1-8501 
- 0-3324 0- 2998 | 02678 
o, | 0-9350 0-8910 0-9249 
on 0-0007 —0-0002 —0-0001 
G3 -0-0005 00005 —0-0004 
264 | 0-3613 | 03613 | 0-3613 
G; 0-0000 00000 | 0.0000 


We come now to the determination of the shapes of the thin lenses. 
Assuming that the two doublets will be cemented, the shapes of the lenses 1, 2 
of each doublet must satisfy the relation 

©1(81—1)(W1—1) = 92(82+1)/(W2—-1)- eee (11) 
There are thus only three independent shapes, say, Sass bs and Six The con- 
dition for zero third-order astigmatism is now set up using equation (44) of 
Part I, the diaphragm being assumed to be coincident with the central negative 
u . 
lens b. In the case of y=—0-375, for example, we obtain 
2.-8746 82, —6 -10688a1 +Se2 +3 44018 ,2+-0-11992=0. 

dition for zero third-order distortion is derived from equation 
(45) of Part I, and for ,= —0-375 we find 

_4-31488%, +-14-18558q1 +Seo+5 84475, —4 8468=0, 


the solution of these two quadratics giving 
§,1=0°22759, S,,=0-29975. 
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Equation (11) then gives S,.=0-42134, and the corresponding relation for the 
back doublet gives S,,——1:92195. The condition for zero third-order coma 
is then written down from equation (43) of Part I and this gives S, directly. 
As a check on the calculations up to this point the third-order aberrations are 
calculated for the thin systems thus obtained. The computations up to this 
stage are summarized in the third and fourth blocks of Table 1. Assuming a 
quadratic relation between the third-order spherical aberration o, and the 
corrector power x, the three values of o, in Table 1 give 


2-4896y2+1-9074y +1 -2562=a,, 


which should be roughly plotted as a guide to the possible correction of spherical 
aberration. 

The next stage involves the replacement of the thin lenses by lenses of finite 
thickness. To do this we introduce axial thicknesses sufficient to permit a 
reasonable aperture, say, f/2-°3, and hence put t,,—0-130, t,,=0-140, and for 
the negative lenses t,.=t,=t,,=—0:03. The replacements are made in the 
standard way which leaves the refraction of the axial paraxial ray the same as 
in the thin system. The third-order aberrations are now recalculated and attain 
the values shown in the first block of Table 2. It is apparent that substantial 
changes have occurred, necessitating readjustment of the shapes. 


The shapes are readjusted as shown in the second block of Table 2 and the 
resulting aberration coefficients are given in the third block of the table. It 
will be seen that now the shapes have been chosen to ensure a small negative 
residual of primary coma and a sufficient value of co, to flatten the tangential 
field. A small positive residual of primary distortion has been left to offset the 
negative secondary and tertiary distortion coefficients. The adjustment of the 
coefficients o, to o; is practically the same in each case, and the only third-order 
quantity that varies with y is the spherical aberration coefficient. Assuming 
as usual, the quadratic relation between y and o,, we find 


1:3293y?41-7971y +0-7608=0. 


This expression is plotted on a ten-times scale in Figure 5. 


Up to this point all the computation has been done on a desk machine. 
Now the complete coefficients of third, fifth, and seventh orders are computed 
on the digital computer, the process requiring about 5 minutes computing 
time per system. This checks all the earlier third-order work and gives the 
higher order coefficients as well. In Figure 5 the variation of the secondary 
coefficients with y is plotted. Each curve is based on three computed points 
and quadratic variation with y is assumed. In addition the dotted curves of 
100, and 10c, are added. Because of the way the shapes are chosen the remaining 


primary coefficients will be practically constant, as instanced by the plot of 
1065. They are therefore not plotted. 


This diagram should be well studied. It summarizes the correction 
possibilities of this type of construction over a wide range of x values. As y 
approaches zero from the negative side the following properties should be noticed. - 
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(1) Spherical aberration.—The third-order coefficient o, is positive and 
increasing, while the fifth-order coefficient uy, is negative and decreases in absolute 


value: values of y can therefore be found for which the spherical aberration can 
be corrected at different apertures. 


: (2) Coma.—The fifth-order circular coma specified by us, wg is positive and 
is decreasing rapidly as y>-0; it is advantageous therefore to proceed as far in 
this direction as possible. The fifth-order elliptical coma (u,, Us, Ys) is positive 


TABLE 2 
DETAILS OF DESIGN EXAMPLE AFTER THE INTRODUCTION OF FINITE AXIAL 
THICKNESSES 

¥: —0-250 —0:375 —0-500 
cory 0-7701 0:7217 0: 7423 
Oz —0-3078 —0- 2932 —0-: 2886 
G3 —0-1104 —0-1204 —0-1322 
204 0-3759 0:3789 0: 3825. 
C5 0-0235 0-0451 0: 0657 

0-0742 0-1384 0-2091 
oye 0-7037 0-5855 0- 4554 
S, —0:3490 —0:3244 —0-3073 
So —2-0990 —2-1363 —2-1760 
Sis 0:3785 0-3951 0-4128 
O1 0-3946 0-2738 0:1945 
Os —0-0231 —0-0247 —0-0207 
63 —0- 0620 —0: 0604 —0-0600 
26,4 0-3969 0:3971 0-3980 
Os 0-0152 0-0123 0-0112 
[4 —5:-008 —7:168 —9-293 
Le 1-548 2-850 4-489 
Us 1-075 1-930 3:013 
Ug —4:-675 —5-595 —6-485 
Us —2-163 —2-576 —2-967 
Us —2 +935 —3-316 —3-731 
[Ly 2-034 1-959 1-939 
lg 1-273 1-218 1-196 
Lg 0-7955 0:777 0-775 
Y-10 —1-076 —1-054 —1-054 
Ua —0-6291 —0:5948 —0-5678 
Uys —0- 0593 —0:0481 —0- 0342 


7 and, as it changes little with y, nothing much can be done about it. Some 


degree of balance can be obtained by offsetting the effects of these positive 
coefficients by adjusting o, to a negative value. This explains why o, was 
adjusted to a small negative value earlier. 

(3) Curvature of field.—The quantities which must be considered are the 
third-order astigmatism o; and the Petzval curvature o,, the fifth-order 
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astigmatism coefficients 4, 4,,, and then the very important oblique spherical 
aberration terms involving wy, w;, 2. Usually it is the values of ay Use (Le 
which set a limit to the extent of the field of a photographic objective. In the 
development of this example we have already made some adjustment of oe 
field by balancing o, and o, to give an approximately flat tangential field as 
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Fig. 5.—The variation of the twelve fifth-order aberration coefficients 
ty tO Uy, with y are shown for the group of type 212 triplets considered 
in the example. The upper dotted line shows the corresponding variation 
of 10o,, the third-order spherical aberration plotted on a 10x scale. The 
lower dotted line shows the variation of 106.2, 6, being the third-order coma. 
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far as third order is concerned. Figure 5 shows that there is a significant decrease 
IM U4, 5, Ue AS x tends towards zero. It will be advantageous therefore to move 
as far as possible in this direction for field correction as well as for the diminution 
of fifth-order cireular coma. 

Examining the spherical aberration coefficients up to the seventh order, 
it becomes clear that, if we select ~=—0-27 approximately, we can achieve 
spherical correction at an aperture of about f/2-3. We can now read off the 
values of all the aberration coefficients of this system and thus are able to assess 
accurately the correction state of this system without any further detailed 
computation. If this is not quite what is wanted, a slight shift in y may give 


LENS SHAPES 


-O*4 -Or2 
CORRECTOR POWER, X 


Fig. 6.—This figure shows the relations between the corrector 

power y and the shapes of the thick lenses ay, b, and c, which give 

simultaneous correction of third-order coma, astigmatism, and 

distortion in the type 212 triplet systems considered in the 

example. The shapes of the lenses for any selected ¥ value may be 
read directly from these graphs. 


a more useful system. The point to emphasize is that Figure 5 and a similar 
diagram showing the corresponding trends among the seventh-order coefficients 
summarize completely the possibilities in this type of system with this selection 
of glasses and residuals. If, for example, the work were repeated with some 
other set of glasses, a corresponding diagram would result, and a direct com- 
parison of the two diagrams would show what advantages or disadvantages had 
resulted from the changes introduced. In this way areas of interest can be 
explored with precise knowledge of what is happening to all the aberration 
coefficients and why the changes are taking place. 
| Suppose now that it is decided that the system with ,=—0-270 is the 
solution appropriate to a particular purpose or requirement. The detailed 
specification of this system can be obtained with very little further computation. 
The initial solution is computed with this value of x, and this gives 

~,=1-2608, ©,=—2°1367, 9,=1:688% 

t,=0-2249, t,=0-1319. 
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Using the values k,/k,;= —0-54083 for the front doublet and k,/k,=—0-44269 
for the back doublet, the powers of the five thin lenses of the system are obtained. 
Figure 6 shows the relations between the corrector power 7 and the shapes of 
the lenses a,, b, and ¢, for which simultaneous correction of third-order coma, 
astigmatism, and distortion is achieved. These curves are based on the data 
contained in the second block of Table 2, i.e. the final shape values used in the 
three computed examples to adjust the third-order coma, astigmatism, and 
distortion. Quadratic relations between the shapes and 7 have been assumed in 
the curve fitting. Hence the shapes required for the lenses a,, b, c, for the 
y= —0-270 system may be read off directly. The values of the shapes of lenses 
a, and ¢c, follow from the cementing conditions, e.g. equation (11). Assuming 
the same axial thicknesses as used already in the example and the shapes deter- 
mined in the manner just described, the radii of the refracting surfaces of each 
lens may be computed and thus the specification completed. Running down the 
intersections of the curves of Figure 5 with the ordinate y= —0-270, the aberra- 
tion coefficient values expected for the system are then read off. These may 
then be confirmed by a 5-minute computation on the digital computer. 
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THE DAMPING OF WATER WAVES BY SURFACE FILMS 
By R. G. VINEs* 
[Manuscript received September 23, 1959] 


Summary 


Measurements have been made of the stilling of small water ripples by surface 
films of cetyl alcohol. The damping is considerable and is somewhat in excess of that 
predicted by existing theories. Imperfections in the ripples induce a directed surface 
drift (surface mass transport) and it is possible that this is indirectly responsible for 
the extra damping. Under natural conditions surface films not only impede the 
formation of small waves but they are also very effective in damping them out. 


I. INTRODUCTION 

The calming action of oil on a water surface has been known from ancient 
times, but the phenomenon was not studied scientifically until the nineteenth 
century (see, for example, Aitken 1883). Lamb (1932) examined the theoretical 
aspects of the damping of water waves and a section of his work was devoted 
to the influence of surface films on capillary waves. More recently similar 
treatments were put forward independently by Levich (1941) and Dorrestein 
(1951). . 

Experimental measurements by Brown (1936) showed a marked increase 
in the damping of small ripples on water when surface films were present. 
However, no quantitative results were reported except for clean water; since 
the wave frequency was high (~300 c/s) damping must have been pronounced, 
and with surface films the experimental errors were probably large. Similar 
measurements have now been carried out at lower frequencies in the presence 
and absence of cetyl alcohol monolayers. The results are in only fair agreement 
with existing theoretical treatments, and the observed differences are well beyond 
the likely experimental error. 


II. EXPERIMENTAL 

The apparatus used is shown in Figure 1; it was similar to the arrangement 
described by Brown (1936) in that stroboscopic illumination of the waves was 
provided by a small vibrating mirror. 

Light from a ribbon projection lamp P passed through a cylindrical lens 
LE, and a narrow horizontal slit S situated close to the lamp. The concave mirror 
- M, (radius 26 cm), rigidly fixed to one arm of a 50 cycle electric tuning fork, 
was used to produce an image of the slit at the point X, and the light was then 
rendered parallel by a large lens L, separated from X by a distance equal to its 
focal length (~15 cm). A plane mirror M,, attached to a travelling telescope 7, 
now deflected the light so that it fell normally upon the water surface in the 


* Division of Physical Chemistry, C.8.1.R.0., Chemical Research Laboratories, Melbourne. 
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enamelled tank B: from here it was reflected back vertically into the telescope. 
The tank B was ~3 cm deep, 30 cm wide, and 35 em long and it was mounted on 
sponge rubber to reduce the effect of stray vibrations. 

Small waves were generated in the tank by the dipper D which just touched 
the water surface. This dipper was a microscope slide ~73 cm across, and it 
was held by a thin wooden rod directly connected to the tuning fork F with a wire 
(Watson 1901). As a result of this arrangement the frequency of the tuning 
fork was reduced to 49-4 c/s. When the fork was in operation ripples were 
obtained on the water which were exactly synchronized with the movement of 
the mirror M,. 

Through the telescope a series of stationary, narrow lines of light could be 
seen. These were images of the slit produced by the curvature of the ripples 
at their troughs and crests. The images from the troughs were above and those 
from the crests below the water surface, at heights which changed progressively 


Fig. 1.—Diagram of equipment used. B, tank; D, dipper; F, 50 c/s 

tuning fork ; L,, cylindrical lens ; L,, lens of focal length f; M,, Ms, 

mirrors ; P, lamp (18A, 6V); S, slit; 7, telescope; W, wire; X, image 
of slit. 


with distance as the amplitude of the ripples decreased (cf. Watson 1901). In 
the present experiments only the images above the water surface—reflections 
from the troughs—were followed. The telescope was focused and one of the 
slit images centred on the cross-wires : a small pointer mounted on a travelling 
stage was then raised or lowered until it, too, was in focus, and its height above 
the water surface was measured. The telescope was now moved along an integral 
number of wave lengths, refocused, and the height of the images above the surface 
determined again with the pointer. In this way it was possible to observe, 
indirectly, the change in amplitude of the ripples at varying distances from the 
dipper, or wave source. 

In some experiments the dipper was dispensed with and ripples were 
generated with a jet of air from a modified fish-tail burner. With care it was 
possible to stabilize the air stream so that waves of frequency exactly twice 
that of the tuning fork were obtained: under these circumstances stationary 


—— 


DAMPING OF WATER WAVES BY SURFACE FILMS 45 


images were again observed through the telescope. It was more difficult to 
produce air-blown ripples at 49-4 c/s ; however, a few rough measurements were 
made at this frequency. 

The wave lengths of the ripples were determined at the same time as the 
other measurements were made ; over the entire distance of travel of the telescope 
(see Fig. 2) there was no detectable change in wave length. 


III. THEORETICAL 
The amplitude of damped straight ripples originating from a long dipper 
is given by 


where uw is the damping factor and @ is the distance from the wave source at 
which the ripples possess an initial amplitude a. When parallel light falls 
vertically on water ripples of this type, the troughs act as cylindrical mirrors 
and series of focal lines are produced above the surface. The focal length f of 
any one of the mirrors is related to the wavelength 4 and amplitude a of the 
ripple, by the equation (Brown 1936) 


FLW Nel ot ital Bs Aa RO CAAT A Cer care ad (2) 
Thus, from (1), 
iO ijietes (0). eee SS (3) 


where f, is the focal length of a ripple of amplitude a). It follows that the damping 
factor p is given by the slope of the linear plot* 


In f=ps-+constant, 2.2... - sane es: (4) 


Theoretical treatments of damping have been concerned with extended 
wave trains, where all ripples have the same initial amplitude a,, which everywhere 
decreases with time ¢ according to the equation 


< being the modulus of decay. 

Although equations (1) and (5) refer to different physical situations, the 
modulus of decay and the damping factor u are simply related to each other in 
terms of the “ group velocity ” V of the wave train (Dorrestein 1951, p. 272), 


that is, 


According to theory (Lamb 1932, p. 624) the modulus of decay +t, when 
- viscosity is the controlling factor in damping, is given by 


egw JOV ME? 1Gary- ondtluuess «hg ture- (7) 


, it is not wecessary to know the actual position of the wave source, 
7 e 


* To obtain a value for 
All measurements were made at distances 


and this was never recorded in the experiments. ane 
greater than ~5cm from the dipper, because of excessive deformation of the surface in its 


immediate vicinity. 
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here y is the kinematic viscosity and k=27/A, 4 again being the wavelength. 
For water at 20 °C (v=0-01 em2?/sec) this gives 


Fyitl ODA icne «yee eee (8) 


When surface films are present, the extensions and contractions of the 
surface associated with the passage of waves are opposed by the variation in 
surface tension so produced. It is this which is mainly responsible for the 
damping of the ripples. Under these conditions (Dorrestein 1951) 


T,=(8/vk2o)?, 
where 
CTSyheed hes”) Jeeps eee Ree ees (9) 


here g is the acceleration due to gravity, 7’ the ‘“ specific”? surface tension 
(surface tension/density), and k, as before, is 27/A. For water the density 
o is unity and Z” is numerically equal to the surface tension T, which is derived 
from the well-known equation for the velocity of propagation v of surface waves 
under the combined influence of gravity and surface tension (Thomson 1871), 


oz (g/l toa T oy: ae ae aes os ee (10) 


This wave velocity v is equal to the product nA where vn is the wave frequency. 


From equations (7) and (9) values of uw corresponding to 7, and +, may be 
obtained in terms of equation (6). These are: 


Dpsenykal VS | 5 dewide okie ged oie 3 (11) 
ee te) iced be fl Pe or se (12) 


The value of the group velocity V is calculated from the relation (Lamb 1932) 


IV. RESULTS 
(a) Clean Water 

The tank B was thoroughly cleaned and filled with water at ~20 °C and 
it was allowed to overflow gently to reduce the concentration of impurities at 
the surface. Typical results are shown in Figure 2 (a), where the logarithm of f, 
the distance of the focal lines above the water surface, is plotted as a function 
of distance w along the tank (the zero is quite arbitrary and x does not represent 
the exact distance from the dipper). The slope of the line drawn is 0-029/cm 
and this must be multiplied by 2-303 to obtain the value of Uobs (0°06; cm). 
The observed wavelength was 0-595 cm (wave frequency 49-4 c/s) and, as may 
be seen in Table 1, the corresponding theoretical value of py, (cf. equation (11), 
with V=40-9 cm/sec) is 0-055/em. From the present measurements the surface 
tension of water, equation (10), is 72-3 dyne/cm. 

Similar measurements by Brown (1936) give 7=72-1 dyne/cm for the water 
used in his experiments. Corresponding to a wave frequency of 300 c/s 
(A=0-172 cm) his value for pop; is 0-45/em : the theoretical value (equation (11) 
with V=77-0 cm/sec) is 0-35/em. Theory requires that 7, the modulus of 
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decay, be proportional to 22, equation (8), and calculation of the proportionality 
constant gives 1-03 from the present results, and 0-97 from Brown’s results. 
Although these figures are reasonably constant in themselves, they are both 
about 20 per cent. lower than the theoretical value of ~1-25 ; thus the damping 
is somewhat greater than predicted. 
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Fig. 2.The damping of small water ripples as a function of distance «x. 
(a) Clean water surface, (b) with high pressure cetyl alcohol film, (c) with 
trace of ethyl stearate. 


Better agreement with theory was occasionally obtained during the present 
experiments, but in general the value of pops was rather teo high. This may 
have been due in part to the presence of slight traces of surface impurities ; 


‘however, in view of the results of subsequent measurements with cetyl alcohol 
24 


films, it is probably a genuine effect. 


(b) Surface Films 
A surface film was introduced by spreading one drop of a solution vs cetyl 
alcohol in a low boiling petroleum ether. The monolayer obtained es 0 pe 
high pressure (~47 dyne/em) and, as is often the case when spreading ta 
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place from solution, this exceeds the equilibrium pressure of a cetyl alcohol film 
in contact with the solid. A typical result is shown in Figure 2 (6); it may 
be seen that, in comparison with Figure 2 (a), the damping is greatly 
increased, for the slope of the line has now risen to 0-18/em (ops =0°41/em when 
[T—25-5 dyne/em). Other measurements were made at different surface 
pressures, and the results of these and the earlier experiments are summarized 
in Table 1. In this table y, is the theoretical damping factor as derived from 
equation (11) (viscosity damping) and p, is the corresponding value obtained 
from equation (12) (damping due to variations in surface tension). 


TABLE 1 
MECHANICALLY GENERATED RIPPLES (WAVE FREQUENCY 49-4 c/s) 
Surface Wave- Slope of Vie u. (per cm) 
Tension length (log, f)/x (cf. 
L r Graph eqn. (13)) — - 
cic) (ex) Uy te | Uy Ue Hobs 
Clean 72°3 0-595 0-029 1-39 0-055 —- | — 0-065 
water 12 1 0-172* 0°19; 1-49 0-35 == 0-45 
Surface 65-5 0-575 0-11f 1-39 0-06 0°17; | 0-23, 0:25, 
films 58-5 0-555 0-12 1-38, | 0-07 0-18; | 0-+25, 0:27, 
42-5 0-500 0-13, 1-37 0-09, | 0-23 | 0-32, 0-31 
25°5 0-425 0-18 1-35 0-15; 0-32; | 0-48 0-41 
j | 


* Measurements by Brown (1936), at wave frequency 300 c/s. ; 
+ This result was obtained at a considerable distance from the dipper : 
damping was less (cf. Section IV (d)). 


nearer to it the 


It may be seen from Table 1 that, when surface films are present, ups is 
distinctly greater than yu, or yp, in every case: in fact ops is roughly equal to the 
sum of yp, and py, over the experimental range. This suggests that damping 
effects due to bulk viscosity and to surface films are largely independent and 
that the total damping is a combination of them both. However, such a view 
receives little support from Dorrestein’s treatment, in which damping by surface 
films is shown to be associated with changes in the nature and extent of hydro- 
dynamic oscillations beneath the surface. Thus the relation u,+p.~ pops 
is probably coincidental. 

According to Dorrestein the most general situation is that in which the film, 
besides showing variation of surface tension with area, also possesses an appreci- 
able surface viscosity ; under these conditions the total damping factor is 
sometimes considerably greater than u,. Nevertheless, calculation shows that 
in the present experiments any such effect is very small indeed,* and it is quite 


* Dorrestein indicates that y. may be greater than py. if N>0-06, where N is a dimensionless 
number related to the surface compressibility. The compressibility of cetyl alcohol films has 
been studied by Nutting and Harkins (1939), and in the present measurements the values of N 
range from about 3 to 40; however, these values are so large that u. exceeds Us by only a small 
amount. If the surface viscosity is also taken into account . is correspondingly reduced, and its 
final value is almost indistinguishable from py. 


tte, > ee 
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clear that the observed damping is always in excess of the theoretical value, 
although the discrepancies are never very large. In Table 1 the difference 
between pp; and yu, is consistently equal to ~0-08; this suggests the presence 
of an extra damping effect which is not included in the theoretical treatment. 


(c) Wind-blown Ripples 

Results for ripples produced by the small air jet are shown in Table 2. 
The surface tension as derived from equation (10) (A=0-535/2 em when 
N=98-8 c/s) is ~28 dyne/cm, corresponding to a film pressure of ~45 dyne/cm. 
This is somewhat lower than the maximum pressure obtained previously, and 
the air blast is possibly responsible ; nevertheless, the film pressure is still high 
and equal to the equilibrium pressure of solid cetyl alcohol. It will be observed 
that uops 18 again greater than y, or uw». However, as in the previous experiments 


TABLE 2 
WIND-GENERATED RIPPLES : CETYL ALCOHOL FILM AT FULL PRESSURE 
Wave Surface Wave- Slope of Vie wu. (per cm) 
Frequency| Tension length (logy) f)/x (cf. 
a. sd X Graph | eqn. (13)) 
(c/s) eyeeiew) (ama) 41 (2 Yate obs 
98-8 ~28 0-535/2 0-32 1-44 0- 0-53 0-82 0-74 
49-4 assumed 0-44* 0:16 1-36 0-1 0-29 0-42 0-37 
28 


* This value may be inaccurate since it was difficult to stabilize the waves in these low fre- 
quency measurements. The wavelength was not determined directly, but calculated on the 
assumption that, as before, 7728 dyne/cm. 


at high surface pressures, pops is rather less than the sum p. +e, the discrepancies 
being roughly the same in both cases. Again, as before, pops—u2=0:08 for the 
experiment at 49-4 c/s, but the difference is much more pronounced at 98 “8 e/s: 
if some extraneous damping effect is responsible it increases markedly with wave 
frequency. 

No measurements were carried out on clean water since there were traces 
of oil in the compressed air line, which quickly contaminated the surface. 


(d) Films of Very Low Pressure 

According to Adam (1941), ethyl stearate gives rise to a gaseous film if the 
surface pressure is less than ~0-03 dyne/em. It was of interest to determine 
whether films of this kind could exert any damping effect, and a minute trace of 
ethyl stearate was therefore introduced to a clean water surface on which ripples 
were generated with the dipper, as before. A small but distinct change was 
immediately observed and the damping was spores: . 20 

The results are shown in Figure 2 (¢c), which, this time, is in the form 
of a curve, ie. the damping is more pronounced as the distance from the 
dipper becomes greater. Similar behaviour, though much less marked, is 
observed with a low pressure cetyl alcohol film (see Table 1). The effect is 


D 
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probably due to movement of the film by surface mass transport, whereby the 
concentration of the monolayer is decreased in the vicinity of the dipper and 
increased towards the far end of the tank. Such movement induced by the 
waves undoubtedly takes place in any film, but its influence is proportionately 
less with films of high surface pressure. 

For the ethyl stearate film the wavelength observed, (A=0-595 cm at a 
wave frequency of 49-4 ¢/s), is the same as that for clean water (cf. Table 1). 
This again gives 772-3 dyne/em so that the surface pressure must be very 
small. Corresponding figures for the damping factor are u,=—0-055/em (from 
equation (11)) and y,=0-16/em (from equation (12)). The mean slope of the 
curve ¢ in Figure 2 is ~0-05 (ups0°12/em), and that of the asymptotic dotted 
line is ~0-07 (uopse0-°16 cm); in either case pytps>Uops> pi It follows 
that the damping effect is quite pronounced, but rather less than is found with 
denser films. This is a reasonable result. 

One interesting experiment was carried out after a film of unknown impurities 
had collected overnight on a clean water surface. The film had a definite 
structure and was quite rigid, and even though its pressure was fairly low 
([=~66 dyne/cm) the damping was pronounced. After some time damping 
decreased and the surface layer could be seen breaking up under the continued 


action of the ripples. This shows that films of high rigidity exert a strong damping 
effect. 


V. CONCLUSIONS 

In this work the amplitude of the individual ripples is small, being of the 
order of 10-% to 10-4 em ;* no results are available for larger waves. Theoretical 
descriptions of the damping of these capillary ripples are confirmed to only a 
limited degree, for the observed damping is always somewhat in excess of the 
predicted values. This suggests that the theories are not completely adequate 
and that some additional damping mechanism is involved. Circular waves are 
damped more rapidly than straight waves ; however, any extra damping due to 
curvature of the wave front must be small in the present experiments. The 
tank used was not large, but the depth of water was always greater than five 
wavelengths so that, for ripples of such small amplitude, any effects of depth are 
scarcely significant. It is possible that damping could be augmented by the 
general circulation in the bulk of the fluid resulting from surface mass transport 
(cf. Section IV (d)). This arises because the waves are imperfect, and the results 
of surface movement of this kind are not considered in any theory. All the 
factors mentioned above lead to an increase in damping, and depend markedly 


on wave frequency in accord with what is observed in Table 2 ; they may all be 
involved in the measurements. 


The experimental facts are at variance with a suggestion of Dorrestein that 
the stilling of small ripples by a surface film is primarily due to the action of 
the film in preventing ripple formation rather than to a direct damping effect. 
He based this view on his result that for capillary ripples of wavelength 


* From equation (2) a~10-4 cm, when f is ~10 cm. 
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0-1-1-7 cm the ratio of the decay moduli T/T, Was never very large, varying 
only from about 1-6 to 4; this implied, he believed, that the influence of surface 
films on the rate of decay of capillary waves was relatively small. However, 
if the damping is considered as a function of distance rather than time (i.e. if 
the factor p. in equation (1) is thought of, and not the factor + in equation (5)) 
the effect of a film is more obvious: the ratio 2/U, 18, in fact, greater than the 
ratio 7,/t, for capillary waves, since the decrease in surface tension which follows 
the introduction of a surface film also produces a corresponding decrease in the 
wave (and group) velocity. The factor (u,-+y.)/u, is greater still, and from the 
present results this would appear to be a more relevant measure of damping 
efficiencies. 

One additional point should be mentioned. Wavelengths are shortened 
in the presence of monolayers, and since the shorter waves are damped more 
easily a further damping is effected indirectly.* It follows that surface films 
not only help to prevent the formation of small ripples as Dorrestein maintains ; 
they are also very effective in damping them out. 
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* This was strikingly demonstrated when a fragment of solid cetyl alcohol was placed on a 


~ clean water surface: the film pressure rose very slowly and the damping of the ripples increased 


progressively as their wavelength decreased. 
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SUNRISE AND ECLIPSE EFFECTS ON THE IONOSPHERE AT 
BRISBANE 


By G. G. BowMAn* 
[Manuscript received November 30, 1959] 


Summary 

A rotating spaced-loop direction-finding system, located at Brisbane, has been 
operated to investigate sunrise effects by using pulsed 3-84 Mc/s transmissions, (a) at 
normal incidence and (b) at oblique incidence. For oblique-incidence recording the 
transmitter was located at Armidale (bearing 202° and distant 355 km from Brisbane). 

Evidence is presented which suggests the formation (at both H- and F,-layer 
levels) of several frontal irregularities, spaced some tens of kilometres apart, extending 
in directions parallel to the sunrise line and travelling, relative to the Earth, with this 
line. These fronts pass overhead at Brisbane approximately half-way between the 
90 km level and ground level sunrise times. 

In the post-sunrise period, F,-layer spreading on the two-hop trace appears, and 
the spreading width increases for about 2 hr, suggesting an F-layer ripple structure, 
with increasing ripple amplitude. The post-sunrise sporadic # occurrence suggests 
frontal irregularities, lying close to the sunrise line direction, soon after sunrise, but 
indicates a swing in direction as time progresses. 

Rotating-loop normal-incidence transmissions were also used to investigate effects 
due to an eclipse of the Sun, at Brisbane, on April 8, 1959. Post-eclipse H- and F,-layer 
frontal irregularities, oriented in directions close to the line representing the end of 
eclipse in the region of Brisbane, suggest a mechanism operating as the eclipse ends, 
which is similar to that operating at sunrise. 


The possibility, that the ripple structure which produces night-time spread-F at 
Brisbane is generated at sunrise, is discussed. 


I. INTRODUCTION 

Sunrise effects have been investigated using a rotating spaced-loop direction- 
finding system located at Brisbane (27-5 °S.; 152-9°EH.). This was capable 
of determining azimuths-of-arrival for signals received away from the vertical, 
and it recorded simultaneously, on separate displays, two pulsed 3-84 Mc/s 
transmissions. One transmitter was located at Armidale (bearing 202° and 
distant 355 km from Brisbane) and the other on the same site as the rotating-loop 
system. Detailed descriptions of this system and the nature of the records from 
it are contained in another paper (Bowman 1960a). Ionograms and fixed- 
frequency h't records have also been used in this analysis. 


II. H-LAYER SUNRISE IRREGULARITIES 
(a) Oblique-incidence Effects 
(i) Introduction.—During the period of rotating-loop recording in 1958, 
oblique-incidence records were taken during 12 sunrise periods in J uly. Within 
five of these periods, additional traces, of a character not found at other times of 
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the day, were recorded. These traces occurred in a Series, spaced from 50 to 
100 km apart. Each trace decreased range at a considerable rate and moved 
towards a range representing the shortest distance from Armidale to Brisbane 
via the # layer. These traces were found to persist, sometimes for as long as 
30min. When the records were good enough, azimuths-of-arrival could be 
read for these signals. The distribution of these azimuths was centred roughly 
on east. This, and the rapid change of range, suggested an association with a 
sunrise irregularity moving with the speed of rotation of the Earth. Plate 1 
Figure 1, shows the nature of the traces recorded. 


VIRTUAL HEIGHT OF REFLECTION (KM) 


eae 25° 272 
ZENITH ANGLE OF INCIDENCE 


Fig. 1—Diagram illustrating method for calculation of virtual 
height of reflection and angle-of-incidence for 3-84 Mc/s oblique- 
incidence transmissions at 0540 on July 31, 1958. 


(ii) Analysis of Sunrise Effect on July 31, 1958.—The first step in this 
analysis was to establish the true ranges of the satellite echoes recorded. To 
do this, a calculation of the group path of the one-hop F,-layer signal was 

_ necessary, as, because of the method of recording, this has to be used as a reference 
level. Use is made of the fact that the virtual height of reflection for the 
equivalent normal-incidence frequency (f=f ; cos 9) is the same as that for the 
oblique-incidence frequency, f’. The oblique-incidence virtual height of 
reflection at 0540 has been calculated in the following manner. Table 1 shows, 
for various angles of incidence 0, (a) normal-incidence virtual aN) read 
from the ionogram at the normal-incidence equivalent frequency (f=f Cos 0) 
and (b) oblique-incidence virtual heights of reflection from the expression, 
177-5 cot 6 (177-5 km is half the distance between Brisbane and Armidale). 
Figure 1 shows how a plot of these variations will give (a) the pea are 
and (b) the virtual height of reflection, at this time. The ionogram is for 
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Brisbane, and, strictly speaking, the virtual height found (406 km) represents 
the virtual height of reflection of 3-84 Mc/s transmissions from stations, one 
midway between Brisbane and Armidale and the other an equal distance on the 
far side of Brisbane (see Fig. 2 (a)). This virtual height gives a group path of 
886 km. 


TABLE 1 
| Normal-incidence | Oblique-incidence 
Normal-incidence Equiv. Virtual Height | Virtual Height 
Freq. for Angle-of-Incidence 6 Or Reta. of Refin. 
(Me/s) | (km) (177-5 cot 0) 
(km) 
3-84 cos 22°=3:-56 418 443 
3-84 cos 23°=3-53 | 410 | 418 
3-84 cos 25°=3-48 400 381 
3-84 cos 27°=3-42 392 348 


Similar calculations have been made for other times within this period, 
and the oblique group paths plotted against time are shown in Figure 2 (ec) 
(using the upper time scale). Since the rapid changes in the F, layer, around 
sunrise, travel round the Earth with a speed equal to that of the sunrise line, the 
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GROUP PATH (KM) 


700 


TIMES FOR (a) 0530 
O610 E.S.T. 


TIMES FOR (b) 0540 0550 O600 06 
10 


Fig. 2.—Calculated group-path variation of oblique-incidence 3-84 Me/s transmissions 
between Armidale and Brisbane on July 31, 1958. A, Armidale ; B, Brisbane ; C, mid point 
between Armidale and Brisbane. 


propagation parameters for the Armidale to Brisbane reflection point can be 
taken as the same as those which pertained 24 min earlier above Brisbane (5 min 
of time being the longitude difference between Armidale and Brisbane). Thus, 
in using Figure 2 (c) the lower time scale gives the correct group path for the 
‘one-hop /’,-layer signal (the o-ray has’ been used for all calculations) 
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Table 2 summarizes the observations made on “ sunrise ” traces of July 31, 
1958. 

Neglecting group retardation and assuming the reflection height to be 
90 km, the data of Table 2 lead to locations, in plan, of the reflection points, 
Shown in Figure 3. The choice of 90 km for the reflection height will be discussed 


TABLE 2 
; Group Path of | Additional Path Sunrise Trace Susan race 
Time F, Trace of Sunrise Trace | Group Path ING 
(km) | (km) (km) , 
0542 887 275 | 1162 080° 
0546 865 225 1090 085° 
| 125 | 990 085° 
0550 823 410 | 1233 — 
350 1173 — 
FAR | 1098 = 
200 1023 = 


later. In Figure 3, positions are plotted relative to a plane which has a constant 
orientation with respect to the Sun. Thus, the positions on it of Brisbane and 
Armidale move with the speed of the sunrise line at 90 km (1500 km/hr). With 
this arrangement, any ionospheric configuration moving with the sunrise line 
Should appear stationary on the plot. No azimuth information is available 
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Fig. 3.—Plan position of reflecting surfaces at the 90 km level, around sunrise 
on July 31, 1958. 


for signals recorded at 0550, so loci of possible reflecting points are plotted in 
this case. The information shown on Figure 3 is consistent with reflection 
taking place from a series of fronts nearly parallel to the sunrise line and moving 
with it. The first front passes overhead at Brisbane at about 0601 Eastern 


Australian Standard Time (H.S.T.). 
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(iii) Analysis of Sunrise Effects on July 13 and 19, 1958.—A method of 
calculation similar to that explained in (ii) has been used for the sunrise traces 
on the oblique-incidence rotating-loop record for July 15, 1958. Figure 4 shows 
the group path plotted against time for these traces ; also for the one-hop F’,- 
layer o-ray trace. Azimuth-of-arrival information is also indicated. The first 


1200 


O-RAY ONE-HOP 
800 F2 TRACE 
a 


GROUP PATH (KM) 


a0o0 
0604 0608 o612 0616 0620 0624 


TIME ‘€.S.T. 


Fig. 4.—Illustration of traces from #-layer sunrise irregularities on 
July 15, 1958. 


sunrise trace terminates at a group path corresponding to great circle reflection 
from the vicinity of the HF layer. At this point it shows as a well-defined trace 
with no change of group path for about 1 min. This is, apparently, the smallest 
range possible for this sunrise trace, and its separation in range from the one-hop 


*1BRISBANE 0610 


pei 
Ca 


ARMIDALE 0610 0620 


Fig. 5.—Plan position of reflecting surfaces at the 90 km level, around sunrise on 
July 15, 1958. 


F,-layer trace (o-ray) allows a determination of the height of reflection of the 
sunrise trace. This was found to be (88--5) km. For the plan-position plots 
(e.g. Fig. 3) the value, 90 km, was assumed. The observed group paths were so 
large that a few kilometres error in height around this level would not affect 
the calculated positions of reflection very much. 


SUNRISE AND ECLIPSE EFFECTS ON THE IONOSPHERE AT BRISBANE 57 


_ Figure 5 is a plot (similar to that of Fig. 3) showing the calculated plan 
positions of reflecting points determined from the ranges and azimuths shown on 
Figure 4. In this case the first sunrise front passes overhead at Brisbane at 
0616 E.S.T. Of the reflection paths shown in Figure 5, 50 per cent. are such as 
would be expected if specular reflection occurred from irregularities lying along 
the sunrise line. i 

Similar calculations for sunrise traces on J uly 13, 1958, indicate the initial 
frontal irregularity passing overhead at 6613 E.S.T., followed by two other 
reflecting surfaces spaced roughly 100 km apart. 


TABLE 3 
Time Transit | Ground | 90 Km 
Ist ‘ Sunrise | Sunrise B-A A-C 
Day Irregularity | B | fal (min) (min) 
ee OR Ee (E.8.T.) 
(E.S.T.) | 
| 
| | 
13.vii.58 | 0613 0641 0552 28 21 
15.vii.58 0616 0640 0551 24 25 
31.vii.58 | 0601 0632 | 0546 Bal 15 


Table 3 indicates the relationship of times, for the first sunrise irregularity 
transit, with the times for ground sunrise and 90 km level sunrise, on the three 
days analysed. 

Although normal-incidence recordings were made at the same time as the 
oblique-incidence records discussed in this section, no traces were seen thereon 
which could be attributed to reflection from these frontal irregularities. 


(b) Post-sunrise Sporadic E Orientation 

On certain days in April and May 1959, the rotating-loop system was 
operated from 0500 to 0800 E.S.T., using the normal-incidence transmitter only. 
Sporadic Z traces were sometimes observed, starting after 90 km sunrise and 
lasting for several minutes. They showed a definite azimuth-of-arrival. Plate 1, 
Figure 2, is an example. It has been shown by the writer (Bowman 19605) 
that sporadic EZ patches are probably frontal in character. Therefore, these 
azimuth readings should give an indication of the orientation of the sporadic # 
fronts. Table 4 shows the days of recording, together with information on these 
sporadic EH occurrences. These short-period sporadic # traces appear to be 
present on most days. 

Baral (1955) has reported, from Calcutta, the occurrence of a ‘‘ sunrise-time ” 
E,, appearing generally between E-layer and ground sunrise times. These 
traces occurred regularly (in more than 60 per cent. of the total observations). 
Figure 6 shows a plot of the time separation between this sporadic # 
~ occurrence and the 90 km level sunrise time, against the angular separation 
~ between the normal to the sporadic # front and the perpendicular to the sunrise 
| line. This reveals that sporadic H fronts occurring soon after the 90 km level 
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sunrise are nearly parallel to the sunrise line. Although the information is 
meagre, the points on Figure 6 representing sporadic H# occurring more than 
1 hr after 90 km level sunrise suggest a swing of this orientation as time progresses. 
This indicates a swing of approximately 90° in 23 hr. 


TABLE 4 
Sunrise Sporadic | Aginntheot Perpendicular | Sporadic E 
Day 90 km H Arival | to Sunrise Duration 
| Level Occurrence | Line (min) 
| 
1959 
April 

8 0526 0538 080° 083° 6 

9 0527 —_ nae = — 
10 0528 —- — — — 
11 0528 0800 180° 082° 10 
12 0529 0546 100° 082° 3 
13 0530 0548 070° 081° 2 
16 0533 0600 080° 080° 5 
y/ 0534 ~- — — _ 
18 0534 0554 090° 080° 7 
19 0535 0642 120° 079° 7 
20 0536 0800 170° 079° 12 

May 

21 0554 0716 100° 070° 10 
22 0555 — as = ee 
23 0556 — — == 
24 0556 0656 090° 070° 3 


III. #,-LAYER SUNRISE IRREGULARITIES 
(a) Pre-sunrise High Multiples 

A regular feature of Brisbane ionograms recorded just before sunrise is 
the presence of high-multiple traces (tenth hop or greater). Plate 2, Figure 2, 
shows such ionograms and Plate 2, Figure 1, an example on an h’t record (at 
2-28 Mc/s). The probability that these traces are due to combinations of H-layer 
and F,-layer reflections is ruled out, since the spacings between the several 
high-multiple traces are approximately equal to the one-hop Ff, virtual height. 
(This is shown most clearly in Figure 7 (a) which represents a tracing of part of 
Plate 2, Figure 1.) The duration of the high multiples which occur at this time 
is seldom greater than 30 min. 

Plate 2, Figure 3, shows a recording of this phenomenon on a rotating-loop 
record indicating an azimuth-of-arrival of 090°. This is consistent with the 
irregularity causing the high multiples being frontal in nature and oriented nearly 
parallel to the sunrise line. The average azimuth-of-arrival for six readings 
taken in April 1959 was 088°, while the perpendicular to the sunrise line during 
this period was at approximately 080°. The average azimuth-of-arrival for 


four readings taken in July 1958 was 082°, and the perpendicular to the sunrise 
line was at 070°. 
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nO ee been made of the occurrence of this phenomenon on 
ee, aie uring 1956. Figure 8 indicates the diurnal variation of 
ewe ae ae = month by month. Also plotted (to the nearest 10 min) 

ae sia ground sunrise. The time of occurrence of the pre-sunrise 
peak in le distribution appears to move with the time of sunrise through the 
year. This pre-sunrise peak is evident from the results of Baird (1954) Fig ure 9 
shows the aggregate of these distributions plotted relative to the fe of saat 
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Fig. 6.—Graph of post-sunrise sporadic H occurrence. Time 

separation from 90 km level sunrise time is plotted against 

angular separation between the normal to the sporadic # front 
and the perpendicular to the sunrise line. 


sunrise. This shows a peak in the distribution approximately 30 min after 
the 250 km level sunrise, and 50min before ground sunrise. The average 
percentage occurrence of 37 per cent. for this peak indicates that, in 1956, this 
phenomenon occurred regularly. The seasonal variation of percentage occurrence 
of this phenomenon in 1956 is shown in Figure 10, when occurrence is considered 
in the period between the time of ground sunrise and an hour prior to this time. 

Referring again to Figure 7 (a), since the pulse repetition frequency used was 
50 c/s, the position of the ground pulse represents 2998 km virtual range when 
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Fig. 7 (a).—Tracing of part of Plate 2, Figure 1, illustrating the nature of high- 
multiple traces in detail. 


Fig. 7 (b).—Comparison between rate of change of group path for trace A from 
(a) and calculated changes for 10-hop and 11-hop transmissions. 


OCCURRENCE (%) 


0300 0400 OS500 0600 
TIME (E.S.T.) 


Fig. 8.—Percentage occurrence of high-multiple traces, around 
sunrise, for the 12 months of 1956. Times of ground sunrise 
(to the nearest 10 min) are indicated. 
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considering high-multiple traces. The traces marked A appear to correspond 
with the same irregularity, since they are in line, whereas the traces marked B 
Suggest a sequence of irregularities, the reflection point jumping from one 
irregularity to the other as time progresses. This multiplicity of traces is a 
common feature of the phenomenon. It is also illustrated by Plate 2, Figures 2 
and 3. 
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Fig. 9.—Percentage occurrence of high-multiple traces (from 
ionograms) for 1956 when plotted relative to ground sunrise. 


It is possible, in some cases, to estimate which multiple hop the trace 
represents and the position of the sunrise F, irregularity relative to the station 
at the time of recording, by the following method. Using the information from 
ionograms, true heights of reflection for a fixed frequency are calculated 


OCCURRENCE (%) 
tv 
ie) ie) 


Fig. 10.—Seasonal distribution (for 1956) of percentage 
occurrence of high-multiple traces for the hour before ground 
sunrise. 


(Schmerling 1958) at 10-min intervals. It is then assumed that the changes in 
height arose from the movement of a fixed ionospheric configuration over the 
Earth with the same velocity as that of the #, sunrise line ; and a true height 
versus distance diagram is constructed on this basis. ‘The results es ane 
July 27, 1958, are shown in Figures 11 (a) and 11 (0) and the group retardation 
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for vertical reflections is shown in Figure 11 (c). The frequency used here was 
2:28 Me/s. ord 
It is now possible to determine, by trial, the location of a scattering irregularity 
on the configuration such that the group path for the multiple-hop echo, at a 
selected time, would fit the observed results. This is not, however, a simple 
geometrical problem, since allowance must be made for group retardation at 
each reflection from the ionosphere. Since this retardation is relatively small 
(in terms of equivalent path), it is sufficiently accurate to assume it the same for 
oblique as for vertical reflection. A guess must be made as to the number of 
hops involved. In the example used (trace A at 0535 on July 27, 1958, Fig. 7 (a)), 
since the total group path for trace A is just less than 12 times the group path 
for one-hop F,, it is clear that the number of reflections involved must be 11 or 10. 
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11 (a)—Proposed ray paths for 1ll-hop F,-layer transmissions for high-multiple 
traces observed at sunrise on July 27, 1958. 


Fig. 11 (6).—Proposed ray paths for 10-hop F,-layer transmissions for high-multiple 
traces observed at sunrise on July 27, 1958. 


Fig. 11 (c).—Variation of group retardation for vertical-incidence 2-28 Me/s transmissions 
at sunrise on July 27, 1958. 


Solutions, based on these alternative assumptions, are shown in Figures 11 (a) 
and 11 (b), the positions of Brisbane, relative to the configuration, at the selected 
times, being shown by dots on the horizontal axis. Another time, several minutes 
later (0547), was then chosen, and the group paths calculated (on the same basis 
a8 before) to the position of the irregularity for 10- and 11-hop transmissions at 
0535. Figure 7 (b) shows the rate of change of group path, calculated in this 
way, for the alternative transmissions. The rate of change of group path for 
11-hop transmission is very close to that observed for trace A, and therefore 
Figure 11 (a) must indicate the approximate position of the principal F, 
irregularity. This appears to pass overhead at Brisbane at approximately 0600. 
This time is 34 min before ground sunrise. 


(b) Post-sunrise Two-hop Spread-F 
Another frequently recurring feature of ionograms at Brisbane is the presence 
of unresolved satellite traces related to the two-hop F, trace: This type of 
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' spreading commences a short time (within 30 min) after ground sunrise. Plate 3, 
Figure 1, illustrates an example of this. Calculations have been made of the 
average spreading widths (in range) of these Spread-2F traces for March 1956, 
April 1956, and July 1958. These average spreading widths are Shown, by 
Figure 12, to grow from a value slightly greater than the transmitter pulse 
width of 10 km to a width in the region of 40 km about 2 hr later. After this, 
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Fig. 12.—Variation of average range-spreading width of two-hop 
spread-F between 0500 and 1200 for March 1956, April 1956, and 
July 1958. 


there is a gradual fall in average spreading width, at least until 1200 (no mone: 
were taken after 1200). The times of commencement of this phenomenon an 
the times of maximum average spreading widths appear to be closely related to 
the time of ground sunrise (Fig. 12). . . 

In view of the apparent orientation of the sunrise F, irregularities along 
directions close to the sunrise line (reported in Section III (a)); and also the 
sunrise line orientation of sporadic E and its subsequent swing in direction 
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(Fig. 6), the azimuths-of-arrival of these satellite echoes appear to be of interest. 
Attempts were made during April and May 1959 to measure the azimuths-of- 
arrival on the assumption that these ionospheric irregularities were of a nature 
similar to those causing spread-F at night (Bowman 1960a). However, the 
signal strength of these echoes was so low that they were not recorded by the 
rotating-loop system. 


IV. Post-ECLIPSE JONOSPHERIC EFFECTS 
During and after the eclipse of the Sun at Brisbane on April 8, 1959, normal- 
incidence rotating-loop records were made. During the eclipse nothing of 
interest concerning the azimuth-of-arrival of signals was recorded. 
However, from 1600, about half-an-hour after the end of the eclipse, two-hop 
F,-layer spreading appeared, the spreading width gradually increasing until 1830 
when it was 150 km, made up of some 10 or more satellite traces. A plot of the 
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Fig. 13.—Variation of range-spreading width of two-hop spread-F' 
at 3-84 Mc/s for 3 hr after the end of the eclipse at Brisbane on 
April 8, 1959. 


variation of spreading width with time is shown in Figure 13. After 1830, 
the range of the main two-hop trace was so great that the extent of the spreading 
could not be determined because of the limited range of recording. It was still 
present at the end of the rotating-loop recording at 2122. The azimuth-of-arrival 
of these satellite signals did not vary by more than 20° in this period, ranging from 
030° to 050°. Plate 3, Figure 2, illustrates the null at 1756. Figure 14 indicates 
that the spread-F fronts, in this case, are aligned roughly in the same direction 
as the line representing the end of the eclipse at Brisbane at 1527 E.S.T. Under 
normal conditions this orientation of the spread-F fronts is rare (Bowman 1960a). 


No nulls were found in the sporadic # traces present from 1530 until 1550 
and from 1616 until 1644, but the trace present from 1740 until 1820 (illustrated 
in Plate 3, Fig. 2) gave an azimuth-of-arrival of approximately 040°. 


tt i rs 
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From 2012 until the end of recording at 2122, high-multiple traces were 
present on the rotating-loop record, indicating echo ranges between 3000 and 
3600 km. These traces showed little change of range during this time, but they 
did show a change in appearance. They first appeared as clean traces, but by 
2122 they showed a number of satellite traces. The azimuth-of-arrival of the 
high-multiple traces was consistent with that found for two-hop spread-F and 
sporadic # in this period (between 030° and 050°). Plate 3, Figure 3, illustrates 
the nature of high-multiple traces at 2120. 


—L 


| 
| 
| 
| 
7 
I 
I 


Fig. 14.—Illustration of times for the end of the Sun’s 

eclipse over Australia on April 8, 1959. A, Eclipse ends 

1400 E.S.T. ; B, eclipse ends 1500 E.S.T.; C, eclipse ends 
1600 E.S.T. © 


VY. DISCUSSION 
The results of this present investigation suggest that, at a time approximately 
half-way between 90km sunrise and ground sunrise, a series of frontal 
irregularities, having spacings of the order of tens of kilometres, are created at 
both the H- and F.-layer levels. These have been detected in the E layer by 
oblique-incidence recording, and in the F, layer by high-multiple traces wee 
appear on normal-incidence records about an hour before ground sunrise. hese 
frontal irregularities (# and F’, layer) are shown to be oriented close to directions 
parallel to the sunrise line and appear to move with it. i 
Evidence is presented to show that, after the occurrence of these ie da 
(moving with the sunrise line) in the # and F’, layers, irregularities of the rons orm 
appear to remain in both layers. Figure 12 shows that, if the rites 
irregularities be regarded as ‘‘ ripples ” in the ionization contours, the was ae 
of these ripples probably increases after sunrise to a maximum ae 4 ieee 
later, since the spreading width appears likely to be a function of ah e amp Ae : 
(The conditions for the occurrence of two-hop spreading, when one- 10p igen : : 
is absent, are discussed by Bowman (1960a).) Evidence for iain it a ne 
irregularities, aligned close to the sunrise line, is obtained from the small p 
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of sporadic H recorded soon after the sunrise period (Fig. 6). Baral (1955) has 
reported isolated sporadic # occurrence, mainly during the sunrise period. 
Bracewell et al. (1951) and Lauter and Schmelovsky (1958) have detected 
sunrise effects on very-low frequency transmissions, which may possibly be 
related to the results of this paper. Recording the signal strength of 16 ke/s 
oblique-incidence transmissions, Lauter and Schmelovsky (1958) found a sudden 
reduction (of the order of 20 dB) in signal strength at the time for 90 km level 
sunrise, at the mid point of the transmission path. A minimum of signal strength 
occurred between 90 km level sunrise and ground sunrise (for the mid point of 
the transmission path). This minimum was followed by a broad peak in signal 
strength (about ground sunrise), after which it assumed a normal day-time value. 


DOWNWARD TILT 


S.R. 250 KM LEVEL S.R. 90 KM LEVEL S.R. GROUND 


Fig. 15.—Variation of F, layer tilt around sunrise at Slough 
(after Bramley 1953). 


It is interesting to note, that, in the case investigated, of pre-sunrise F’,-layer 
high multiples, the F,-layer sunrise irregularity appears to pass overhead at 
approximately the same time (0600 on July 27, 1958) as the leading H-layer 
sunrise irregularity four days later (0601 on July 31, 1958). This suggests an 
effect generated at the H-layer level and transmitted to the F,-layer by electro- 
dynamic coupling. 

Figure 15 reproduces a result due to Bramley (1953) showing the variation 
of the F,-layer tilt around the sunrise period. Times for sunrise at 250 km, 
90 km, and ground levels are indicated. Figure 15 reveals a tilt of the F, layer 
which increases with time and commences at the 90 km level sunrise rather than 
at the 250 km level sunrise. This may be important in view of the coupling 
suggested earlier. As pointed out previously, Figure 12 suggests ripples of 
increasing amplitude after sunrise. Therefore the marked oscillation of the F,- 


layer angle of tilt after the initial sunrise effect (Fig. 15) appears to be consistent 
with the results shown in Figure 12. 


Conditions in the ionosphere after the eclipse on April 8, 1959 are very 
similar to those found after sunrise. If the uncovering of the Sun after the 
eclipse be regarded as a virtual sunrise, the line along which the eclipse ends is 
analogous with the sunrise line. The gradual increase in two-hop F, spreading 
widths after the eclipse, and the continued orientation of E- and F,-layer 
irregularities along directions close to the line representing the end of the eclipse, 
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- all point to a mechanism operating as the eclipse ceases that is comparable with 


that operating at sunrise. 

The possibility arises that ripples, generated either at sunrise or sunset, are 
responsible for the one-hop spread-F observed at night at Brisbane. Effects, 
similar to those discussed here for sunrise, have not been detected at sunset, but 
there is some evidence to suggest that ripples, generated at sunrise, persist 
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Fig. 16.—Diagram illustrating comparison between seasonal 

variation of spread-F’ at Brisbane (1944-1958 inclusive) and the 

angular separation between the sunrise line and the magnetic 
meridian. 


throughout the day, and are responsible for the spread-# observed at night. 
The two-hop, post-sunrise spreading width, illustrated by Figure 12, should be a 
measure of the amplitude of the ripples. At 1200, it is obvious that a ripple 
structure is still present, particularly for July, a month of maximum spread-F 
occurrence. Another possible association between ground sunrise and night-time 
spread-F is shown in Figure 16. Here, the seasonal distribution of spread-F' 


Fig. 17.-Seasonal distribution of sunrise high multiples 
for 1956, after an allowance is made for a possible sunspot- 
cycle variation. 


(for the years 1944-1958 inclusive) is compared with the seasonal variation in the 
angular separation between the sunrise line and the magnetic meridian. The 
similarity is interesting, particularly as the spread-F minima are displaced 
somewhat from the equinoxes and occur roughly at times during the year when 
the sunrise line lies along the magnetic meridian. 

Another apparent association between sunrise and night-time phenomena 


_ is revealed by considering the seasonal distribution of the sunrise high multiples. 
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During 1956 the sunspot activity rose sharply, and therefore the distribution 
shown in Figure 10 suggests a direct relationship with sunspot activity. If this 
association is assumed, a seasonal distribution, with some independence from 
sunspot activity, is obtained by subtracting from the percentage occurrence of 
each month, the background variation proposed on Figure 10. Figure 17 shows 
the distribution deduced in this way, revealing equinoctial maxima, a winter 
minimum, and a summer sub-minimum. This distribution is the inverse of that 
found for spread-F' at Brisbane. One possible clue to this inverse relationship 
is the need for a small ripple amplitude to produce high multiples (Baird 1954) 
and much larger amplitudes for one-hop spread-F. 
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A RELATIONSHIP BETWEEN SPREAD-F AND THE HEIGHT OF 
THE F, IONOSPHERIC LAYER 


By G. G. BowMan* 
[Manuscript received November 30, 1959] 


Summary 


The variation, throughout the night, of range-spreading spread-F widths, at 
Brisbane, is investigated. Evidence is presented to support the hypothesis that the 
ripple amplitude of the spread-F ionospheric irregularities varies directly with the 
layer height. This leads to an association between the range-spreading width and the 
height of the layer. 

Certain aspects of frequency-spreading are also discussed. 


I. INTRODUCTION 

In a previous paper on spread-F at Brisbane (Bowman 1960,7 subsequently 
referred to as I) the diurnal variation of spread-F was considered statistically, 
the existence or not of spread-/ traces being taken as a measure of its occurrence. 
No account was made of the widths (in range) of the spread-F traces. Treated 
in this way, the diurnal variations suggested increased spreading occurrence 
when the F, layer increased in height. It was postulated that, possibly, this 
resulted from an increased amplitude for the spread-F ‘“‘ ripples ”’, as the height 


WA 


ee 


Fig. 1—Diagram illustrating suggested dependence of range- 
spreading width on the amplitude of the irregularity ripples. 


of the layer became greater. In this case, it would be expected that the width 
of range spreading, at any particular time, should also be related to the height 
of the layer, for this spreading width is determined by the range of zenith angles 
from which satellite signals are possible, and this range should increase with 
ripple amplitude, as indicated in. Figure 1. The present paper deals with this 


aspect of the phenomenon. 


* Physics Department, University of Queensland, Brisbane. 
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Fig. 2.—East-west profile of height changes and plan position of F, 

irregularities, (a) leading to a maximum range-spreading width some 

time later than time of maximum layer height, (0) leading to a maximum 
range-spreading width at the time of maximum layer height. 
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Fig. 3.—A comparison between h’F, variation and range-spreading 
spread-F' widths, (a) for the night of August 16/17, 1956, and 
(6) for the night of July 23/24, 1956. 
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II. ANALYSIS AND Discussion 

If the hypothesis discussed in Section I is correct, it would be expected that 
a relationship could be established between height changes and the relative 
widths of the range-spreading trace. However, this will be complicated by the 
fact that, in general, the frontal orientation of the ripples is not the same direction 
as that of the gradient of height changes. In I it was established that the most 
favoured azimuth-of-arrival for spread-F echoes occurred in the north-west 
quadrant. Therefore, if it is assumed that height changes have an east-west 
gradient, and the orientation of the ripples is at an angle of 45° to this direction, 
the arrangement is as illustrated in Figure 2 (a). In this figure the height profile 
is imagined as remaining stationary, while the position of the station moves. 
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Fig. 4.—A comparison between f,/, variation, h’P, variation, 
and range-spreading spread-F' widths for the night of May 3/4, 
1956. 


The frontal irregularities are closely spaced, and, as an approximation, they are 
considered to have infinite length. This figure illustrates that, since the ripples, 


which are responsible for the extreme edge of the range spreading, are appreciably 


displaced from the zenith, they will correspond to the layer height over the 
station some minutes earlier. Thus maximum spreading width can be expected 
some minutes after the maximum layer height. This can usually be observed a 
the records. Figures 3 (a) and 3 (b) show examples. For the F,-layer height 
ntity Amin is used. 
A a eae the spread-F' phenomenon at Brisbane for the night of May 3/4, 
1956 is such as to suggest that some of the complicating parameters are ce 
The layer height does not oscillate throughout the night, as it usually does, bu 


- shows a more or less gradual rise to midnight, followed by a gradual fall. The 


peak in range-spreading width coincides with maximum height of the layer, so 
co) 
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possibly the frontal orientation of the ripples, at this time, lies in the same 
direction as the height change gradient. Figure 2(b) shows the suggested 
configuration. Figure 4 illustrates the f,F,, h’F,, and range-spreading width 
variations on this night, the readings being taken from the ionograms made at 
10-min intervals. The range-spreading width variation follows closely the h’F, 
variation. 

It was explained in I that, possibly, when the ripple amplitudes become 
relatively small, the ripples can still be detected in the form of frequency 
spreading, provided there is a sufficiently large f,F, gradient across the ripple 
structure (cf. I for definitions of range spreading and frequency spreading). 
The ionograms on this night are particularly clear in showing the transition 
from range spreading to frequency spreading. It is seen that, by 0300 hr, 
the range-spreading width has been reduced to a low value, but after this time 
the ionograms show increasing frequency spreading. Also, after 0300 the fof; 
value changes at a fast rate. This appears to be evidence for the mechanism, 
proposed in I, to explain frequency spreading. 

Sample ionograms at various times throughout the night of May 3/4, 1956 
are shown in Plate 1. 
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THERMAL TRANSIENTS ASSOCIATED WITH NATURAL CONVECTION 
By R. C. L. Boswortu* and C. M. GRropENt 
[Manuscript received November 6, 1959] 


Summary 


An electrical circuit, equivalent to the thermal processes concerned in the act 
of setting up a system of natural convection in a fluid, has been proposed. This circuit 
consists of parallel elements, one attributed to convection alone and one to conduction 
alone. The former consists of a resistance and an inductance in series, the latter of 
a resistance and a shunt capacity. The nature of the transients associated with each 
circuit has been derived by Heaviside analysis. Functions, at the most involving 
no more than two independent variables, are found connecting easily measured character- 
istics of the transients and the circuit parameters. These functions are, it is proposed, 
to be used to examine experimental thermal transients and through them to examine 
their circuit parameters. 


I. INTRODUCTION 

The way in which the distribution of temperature settles down after a 
sudden thermal shock is ordinarily described by the Fourier equation. In the 
Fourier equation a single characteristic property of the system is referred to the 
thermal diffusivity, which in turn is equal to the ratio of the measures of two 
properties of the medium in which heat flow occurs, namely, the conductivity 
and the capacity per unit volume. For a system of given geometrical dimensions 
the diffusivity gives rise to a single time constant, namely, the ratio of the thermal 
capacitance to the thermal conductance. 

The representation described by the Fourier equation can accordingly be 
expressed, at least in principle, by an equivalent electrical circuit involving only 
resistances and capacitances. Systems which involve convective as well as 
conductive flow may be more complicated. Frank-Kamenetskii (1938), for 
example, claimed that two time constants are associated with thermal convection. 
The nature of the initial heat flow associated with convection, especially in liquid 
systems, is formally quite different from that from a resistance-capacity circuit, 
and one of the authors (Bosworth 1946, 1948) has suggested that these systems 


exhibit a phenomenon properly called thermal inductance. These phenomena 


are connected with the temperature changes which follow the sudden steady 
application of heat to a wire immersed in the fluid. The heat flow from the wire 
to the fluid was constant, while the temperature difference between the wire and 
the fluid, which was initially zero, finally settled down to a steady value after a 
sufficient time of heating. In the phenomenon that one of the authors described 
as thermal inductance the temperature difference at first rose to a value higher 


than the final value and approached it from above. 


* School of Chemistry, University of New South Wales, Sydney. 
+ School of Mathematics, University of New South Wales, Sydney. 
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In certain fluids, notably in gases and in very viscous liquids, the initial rise 
in the temperature difference above the final value was not realized ; in other 
words, the final temperature value was approached from below. Such a transient 
may be regarded as of a purely capacitive type. On the other hand the transient 
which approaches its full value from above could be described as an inductive 
type. Thermal transients of both types have been described by Ostroumoff 
(1956) and interpreted in a purely mechanical manner. A more exhaustive 
study has shown that moderately viscous liquids appear to give a capacitive-type 
transient at a low heating rate and an inductive-type transient at a higher heating 
rate. 

Typical transients for each type are shown by a hot wire of 30 S.W.G. 
immersed in an 81-2 percent. (supersaturated) sucrose solution at 18-0 °C. 
Figure 1 shows the thermal transient for a heating current of 0:52 A and Figure 2 
shows that for a heating current of 2:65 A. The former appears to be of the 
capacitive type, the latter of the inductive type. 


II. AN ASSUMED EQUIVALENT CIRCUIT 

On the application of the flow of heat from the wire to the fluid an initial 
effect results in setting up a temperature difference between the wire and the 
outer regions of the fluid, and later also a temperature difference occurs between 
the fluid in close contact with the wire and the more distant regions. At the 
same time, there exists initially no circulating system of convection currents 
and the thermal resistance opposing heat flow is initially high. When the system 
of convection currents has been fully developed the resistance to thermal flow is 
considerably reduced. The equivalent circuit representing the heat flow thus 
must involve two resistances in parallel, one of which r is always operative and 
the other R only becomes operative after the system of convection currents 
has been developed. In the equivalent circuit we represent # as in series with 
an inductance L. The whole system has an additional capacitance C’ which 
operates in parallel to r and R+L. This takes into account both the heating up 
processes and the action of setting up a steady temperature gradient through 
the fluid. The capacitance C is expected to include the thermal capacity of the 
wire and of the fluid adjacent to it. 

The overall suggested equivalent circuit is represented in Figure 3. The 
circuit, as an electrical one, has been discussed by Carter (1944) in connexion 
with the transient voltage across the Tirrill regulator contacts on a generator. 
Carter expressed the input impedance in terms of Heaviside operators and 
interpreted the operators in terms of exponential functions. For certain ranges 
of the variables (r, R, 0, and L) the exponents become complex quantities and the 
transient takes the form of a damped oscillation. When the exponents are both 
real the disturbance is overdamped and there are no oscillations. 


By matching actual values of the thermal transient to the theoretical values 
of the calculated transient with selected values of &, C, L, and 7, it is, at least in 
principle, possible to derive all these associated cireuit characteristics. In 

practice it is relatively easy to match the circuit characteristics from the oscillatory 
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transient. Points may be selected from values of the temperature which pass 
through the asymptotic value and which attain one or more turning points. 
However, the derivation of such characteristics as these from the theoretical 
expressions for the circuit diagram is not particularly convenient in the form given 
by Carter, and an equivalent, but more convenient, expression may be derived. 


Fig. 3 


III. THE CHARACTERISTICS OF THE HQUIVALENT CIRCUIT 
The overall formal impedance Z of Figure 3 may readily be seen to be 


= r(R+Lp) 
OLrp?4-ROri4- Lp 47 Lie alee (1) 


where p is the Heaviside operator 


Let the current take the form of the Heaviside function, namely 


I=0, when t<0, 
I=1, when t>0. 


The potential for any positive value of ¢ then takes the form ZI » where 


peli) ‘1 ee (? Ror) : oosn V{(L—ROr)?—4°0L}, 


rR 2rCL 2rCL 
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This equation may be written in a more concise form by means of the 
following substitutions : 


X= (G RO) 2raf (DO), a Sos ors cov es (2a) 

Vet LC DAT ML Ge) one Sirti ee teehee ee (2b) 
and 

eet GA) MS beeen cease oes oe (2c) 


Equation (2) then becomes 


_ rk 
r+R} 


Z 1—exp—(X+y}s| cosh 4/ (X? —1)t— 


Vi xa 1) 


When the modulus of X is less than unity the arguments of the hyperbolic 
function become complex and may be replaced by trigonometrical functions, 


meth, ( 1 X+Y 
= Rt —exp—(X+7} | cos a/ (1 —X?)t — V(t —X3) sin 4/(1 —Xij| t 


valid when —1<X<1. 


The oscillatory type of transient is represented by equation (4). A com- 
pletely overdamped transient is represented by equation (3) with the value X 
negative and numerically greater than unity. The positive values of X that 
are greater than unity show no true oscillatory motion. However, the expression 
in the square brackets in equation (3) can always change from a positive to a 
negative value with an increasing +; this means that the transient value of Z 
overshoots the asymptotic value and finally approaches that value from above 
and not from below as in the region for which —1<X<-+1. 

Experimentally we may expect a capacitive type of transient to appear 
when X<—1 and an inductive-type transient when X>—1. In the limited 
region —1>X>1 the transient may show a series of damped waves. However, 
in all thermal flow circuits we expect the damping to be fairly severe and any 
wave other than the first may well be swamped by the noise level. 

As indicated earlier the four most readily attainable characteristics of 
experimental ‘“ inductive-type ” transients are the value ¢, of the time required 

to first attain a value equal to the asymptotic value, and the value ¢, to attain 
the first stationary point and the corresponding impedances Z, the asymptotic 
value, and Z,, the value at the first maximum. Since four variables are necessary 
to describe the transient, the four easy experimental values are, at least in 
principle, sufficient to fix the nature of the transient. 


The value of ¢, from equation (3) takes the value 


n 1 /(X?—1) 
Ye 


/ (LC) (KX? —1) are tanh 


=a) ATG CORD Xief) Geis. 6 o aye poms « (5a) 
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From equation (4) the corresponding value reads 


eae are con Al?, Ges een (6) 


VJ(Le) VW(1—X?*) 


The value of ¢, from equation (3) takes the value 


tee V(X?—1) 
VULO) ~ V(k?=1) are tanh Tava Pee (7) 
and from equation (4) 
fo zZ arate! eee (8) 


/(0) 4/A—X?) "| yeh os Se hs wr aLetane 


The ratio of these two times, which we shall take as equal to A, becomes, 
in the trigonometrical region, 


t, are cos X 


i) are tan (V—X)(X+Y)} > 


We now make a substitution for X, putting 


XA=c08 T; 
and obtain 
_sin (1—1/A)T 
atid (9) 
In the hyperbolic region a corresponding substitution, 
X=cosh 7, 
gives 
_ sinh (1—1/A)T 
emt nial! Ta 2 oe ee (9a) 
Let the ratio of the corresponding two impedances be 
Z,/Z =. 
The trigonometrical region then gives 
X+1/Y 
B=1+Y ex a ee xi, ee ats 1 
and, in the hyperbolic region, 
X+1/Y ; 
B=1+Y ex ie ecce 2a are cosh ae ripe 1 
EYE = et 
Substituting T for X we get 
qr i 
In Yi Bj] = 
n Y/(B—1) in (<0 T+3) Bahama vos (11) 
or 
In Vie cosh p+ Li 
ainh 7 4 (11a) 


in the respective region. 
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» Elimination of Y between equation (9) and (11) now gives 


cee 
In (B—1)=In ae ma | Tecot- (1 AV hae ne (12) 
or ; 
inh (1—1/A)T 
in(B—1)=In aaa ra | T coth (1—1/A)T, .. (12a) 


in which one of the required combinations of the circuit characteristics (namely, 
T=arc cos X or are cosh X) are implicitly expressed in terms of measurable 
characteristics of the transients, namely, A and B. 


IV. EVALUATION OF CIRCUIT PARAMETERS FROM THE CHARACTERISTICS 
OF THE TRANSIENT 

The relationship between X and A and B can be obtained from tables, 
such as Table 1, derived from equations (12) and (12a). Reading off from 
values of A and B (interpolated, if necessary) it is possible to derive X and thus 
(L—RCr)/2r,/(LC) from Table 1. 

Table 1 in effect gives a relationship between the parameters A and B at a 
series of constant values of XY. A more convenient tool is a relationship between 
A and X at a series of fixed values of B. Such data, derived from Table 1, are 
shown in Figure 4. It will be noted that the curves showing A versus AX at 
constant B are practically rectilinear at the higher values of XY. For limiting 
high values beyond the range covered by Table 1 or Figure 4, equation (12a) 
may be written in the approximate form 


e(—1/A)T 
In (B—1)=n —qqz —T; 


or 
Bae 7ti4 
=a( 2X) Ae rete rien en ane (13) 
Under similar conditions the limiting value of Y becomes 
GD. ae) re en (13a) 


so that the apparently linear relationship between A and X does not apply 
to still higher ranges. 


V. Tue Ratio oF CONDUCTIVE TO CONVECTIVE RESISTANCE 

By using Figure 4 and equations (9) (or (9a)) it is possible to derive both 

X and Y from the known values of A and B. Now X and Y are connected 
from their primary definitions in equations (2a) and (2b) and, in fact, 


OXY ALY Ry. vs eon et es (14) 
Accordingly, the ratio R/r thus depends only on the parameters A and B. 
The relationship proves to be a particularly simple one. Substituting equation 
(9) for equation (14) gives 
2 cos T sin T/A sin? T/A af, 
“sin (T—T/A) ‘sin? (L—T/A) ’ 
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’ TaBLe 1 


VALUES OF B FOR GIVEN X AND A, COVERING BOTH TRIGONOMETRICAL AND HYPERBOLICAL REGIONS 


- 
\ A 1-2 1-3 1-4 1-5 1-8 = 2:0 255 7) 33-0 3-5 4-0 4°5 
ae 
—1:0 1:004+ 1-029 1-082 1-163 1-575 2-000 3-776 7-135 13-248 24-141 43-275 
—0:9 1-002+ 1-015 1-046 1-093 1-324 1-539 2-276 3-297 4-619 6-258 8-237 
—0:8 1-002— 1-012 1-037 1-077 1-264 1-435 1-991 2-706 3-567 4-563 5-687 
Nor 1-QO0L-- 2-011 1-032 1:066 1-228 1-373 1-831 2-397 3-051 3-779 4-573 
06 1:001+ 1-0095 1-028 1-059 1-203 1-330 1-726 2-199 2-782 3-312 3-933 
{1005} 1:001+ 1-009 1-026 1-054 1-184 1-298 1-648 2-059 2-513 3-000 3-513 
—0-4 1-001 1-008 1-024 1-049 1-169 1-273 1-588 1-953 2-351 2-773 3-213 
—0:3 1-001— 1-007 1-022 1-046 1-156 1-252 1-540 1-870 2-226 2-600 2-987 
—0-2 1-001— 1-007 1-021 1-043 1-146 1-235 1-501 1-803 2-126 2-463 2-811 
—0°1 1-001— 1-006 1-0195 1-040 1-137 1-221 1-468 1-747 2-044 2-352 2-669 
0 1-001— 1-006 1-018 1-038 1-129 1-208 1-440 1-699 1-975 2-260 2-551 
0-1 1-001— 1-006 1-017 1-036 1-122 1-197 1-415 1-659 1-916 2-181 2-452 
0-2 1-001— 1-005 1-017 1-034 1-116 1-187 1-393 1-623 1-865 2-114 2-368 
0:3 1-001— 1-005 1-016 1-033 1-111 1-178 1-374 1-592 1-821 2-056 2-294 
0-4 1-001— 1-005 1-015 1-031 . 1-106 1-170 1-357 1-564 1-781 2-004 2-231 
0:5 1:001— 1-005 1-0145 1-030 1-101 1-163 1-342 1-539 1-746 1-959 2-175 
0-6 1:001— 1-005 1-014 1-029 1-097 1-157 1-328 1-517 1-715 1-918 2-124 
0-7 1-:001— 1-004 1-013 1-028 1-094 1-151 1-315 1-497 1-687 1-882 2-078 
0-8 1-001— 1-004 1-013" 1-027 1-090 1-145 1-304 1-479. 1-661 1-848 2-038 
0-9 1-:001— 1-004 1-0125 1-026 1-087 1-140 1-293 1-462 1-638 1-817 2-002 
iol 1-0005 1-004 1-012 1-025 1-085 1-135 1-282 1-446 1-616 1-792 1-968 
1-2 1-0005— 1-0035 1-011 1-023 1-079 1-127 1-266 1-418 1-579 1-742 1-909 
1-4 1-0004+ 1-0035 1-011 1-022 1-074 1-120 1-251 1-395 1-546 1-700 1-857 
1-6 1-0004+ 1-003 1-010 1-021 1-070 1-113 1-237 1-374 1-518 1-664 1-814 
ihote 1-0004— 1-003 1-0095 1-020 1-067 1-107 1-225 1-356 1-493 1-633 1-775 
2-0 1-0004— 1-003 1-009 1-019 .1:064 1-102 1-215 1-339 1-470 1-605 1-741 
2-5 1-0003+ 1-0026 1-0081 1-017 1-057 1-091 1-193 1-305 1-424 1-547 1-672 
3:0 1-0003+ 1-0024 1-0073 1-015+ 1-051 1-083 1-175 1-279 1-388 1-502 1-617 
3:5 1-0003— 1-0022 1-0067 1-014 1-047 1-076 1-161 1-257 1-359 1-465 1-573 
4-0 1-0003— 1-0020 1-0061 1-012+4 1-043 1-070 1-149 1-239 1-335 1-435 1-537 
4-5 1-0002-+ 1-0019 1-0057 1-012— 1-040 1-065 1-139 1-224 1-314 1-409 1-506 
5-0 1-0002+ 1-0017 1-0053 1-011— 1-037 1-060 1-130 1-211 1-297 1 387 1-480 
5°5 1-0002+ 1-0016 1-0049 1-010+ 1-035 1-057 1-123 1-199 1-282 1-368 1-458 
6-0 1-0002— 1-0015 1-0046 1-010— 1-033 1-053 1-116 1-189 1-268 1-352 1-438 
6-5 1-0002— 1-0014 1-0043 1-009— 1-031 1-050 1-110 1-180 1-256 1-337 1-420 
7:0 1-0002— 1-0014 1-0041 1-008+ 1-029 1-048 1-105 1-172 1-246 1-324 1-405 
Tots) 1-0002— 1-0013 1-0039 1-008+ 1-028 1-045 1-100 1-165 1-236 1-312 1-390 
8-0 1-0002— 1-0012 1-0037 1-008— 1-026 1-043 1-096 1-159 1-228 1-301 1:377 
8-5 1-0001+ 1-0012 1-0035 1-007+ 1-025 1-041 1-092 1-153 1-220 1-291 1-366 
9:0 1-0001+ 1-0011 1-0034 1-007— 1-024 1-040 1-089 1-147 1-213 1-282 1-355 
9-5 1-0001+ 1-0011 1-0032 1-007— 1-023 1-038 1-085 1-142 1-206 1-274 1-345 
10-0 1-0001+ 1-0010 1-0031 1-006 1-022 1-037 1-082 1-138 1:200 1-267 1-336 


eS ee ee eee ee ee 
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- with n the ratio of the resistances. This expression may readily be reduced to 


sin? 7 =(1—n) sin? T(1—1/A), ............ (15) 
with a corresponding expression in the hyperbolic region. 


The quantity » may take values from 0 to co but may never be negative. 
The upper limit corresponds to a type of heat flow in which the contribution 
by the convective transfer has vanished. The lower limit corresponds to a 
condition at which the convective contribution is infinitely large in comparison 
with the conductive transfer. The only solution of equation (15) pertinent to 
practical problems for which n vanishes is the particularly simple form 


ASTiOT—x),, 4m Bn beeen: (16) 


The points where the curve of R/r—0 on Figure * cross the Lees of constant B 
are solutions of equation (16). The region to the aene of this line corresponds 
to finite values of R/r, that to the left to unattainable values of R/r. Any 
experimental value of A and B lying to the left of vie line thus porespon eas a 
condition which is inconsistent with the equivalent circuit postulated in Figure 3. 

Values for A at a range of values Ob F (Or) may be obtained from equation 
(15) for selected values of n (0-1, 0:2, etc.). The lines at selected values of 3 
are drawn a8 a series of contours crossing the constant B lines in Figure 4. 


- Accordingly, from this figure we can read off, given a known pair of values for 
A and B, the corresponding values for X and for the ratio R/r. 
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VI. TREATMENT OF EXPERIMENTAL RESULTS 

Experimental results are obtained in the form of a transient curve plotting 
the temperature excess (or alternatively the quotient of the temperature excess 
by the heat flow rate) against the time. The peak temperature and the time to 
attain that peak are recorded, together with the final temperature and the time 
first to cross that amplitude. From these measurements the ratio of temperatures 
(B) and the ratios of the times (A) can be read off. From Figure 4 we may now 
read off the value X and also the ratio R/r. Alternatively, from the X value 
we may compute the value of Y using equation (9) or (9a)). From Y and X 
the value of R/r can also then be derived by use of equation (14). From the 
final temperature 

Zi Erf(R--7), sss doassewee eon (17) 


it is then possible to obtain values for both R and r. 


TABLE 2 
VALUES OF THE RATIO 4/(LC)/t, AS FUNCTION OF X 


x Vv (LC)/ty x Vv (LC)/ty x Vv (LC)/t 
—1:0 0-0000 0-6 0-8627 5-5 2+263 
—0:9 0-1620 0-7 0-8978 6-0 2-387 
—0:8 0+ 2402 0-8 0-9324 6-5 2-510 
—0-7 0-3044 0-9 0-9665 7-0 2-630 
—0:6 0-3613 1-0 1-0000 7-5 2-749 
—0:5 0-4135 1-2 1-066 8-0 2-866 
—0:4 0-4624 1-4 1-130 8-5 2-983 
—0°3 0-5086 1-6 1-193 9-0 3-098 
—0:2 0-5529 1:8 1-255 9-5 3-211 
—0-1 0-5954 2-0 1-315 10-0 3-324 

0 0-6366 2-5 1-462 12-0 3-763 

0-1 0-6766 3-0 1-605 14-0 4-191 

0:2 0-7155 3-5 1-743 16-0 4-608 

0:3 0+ 7534 4-0 1-877 18-0 5-018 

0-4 0: 7906 4-5 2-008 20-0 5°415 

0-5 0-8270 5:0 2-137 


Again from equation (2b) it is possible to obtain a direct measurement of 
V(L/C), 
/ (L/0)=RY. 


A complementary expression may be obtained for the quantity ./(LO). The 
ratio t,/+/(LC) using equations (5a) and (6), may be expressed as a function of XY. 
Numerical data for this function (covering the range —1<X<10) are given 
in Table 2. 

From this table the values of 1/(ZC) may be derived from 1, knowing the 
value of X. With a very high value of X the expression takes a limiting form 


/ (LC) =Xt,/2 +303 log 2X. 


———— 


oe ee ep re 
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From known values of 4/(L/C) and 4/ (LC), independent values of Z and 0 may 
be obtained, and the identification of all experimental parameters is now 


complete. Examples illustrating these methods will be given in the following 
paper (Bosworth 1960). 
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ATTEMPTS TO MEASURE THE INDUCTIVE ELEMENT ASSOCIATED 
WITH THE NATURAL CONVECTION OF HEAT 


By R. C. L. BoswortTH* 
[Manuscript received November 6, 1959] 


Summary 


A study has been made of the variation in time of the temperature of a wire immersed 
in a fluid and heated by a constant electric current. For a given fluid the curve obtained 
by plotting the ratio of the temperature of the wire to the heat input versus the time 
is initially the same shape for all rates of heat input. Divergences from the lowest 
heating rate set in only when the system of convection currents sets in. This occurs 
at earlier times after the commencement of heating the higher the heating rate. Expres- 
sions already developed are used to evaluate the resistive, capacitive, and inductive 
elements required to fit the observed transient curves. The values of the former two 
types of element are consistent with an assumed stagnant film of a thickness the order 
of 1mm around the heated wire, but the value of the deduced inductive element is 
some 104-105 greater than that associated with the kinetic energy belonging to the 
system of convection currents. 


I. INTRODUCTION 

The work described in the previous paper (Bosworth and Groden 1960) 
has provided a means for evaluating certain thermal properties of a fluid from 
the study of convective heat transfer in the transient state. The object of the 
present paper is to describe some preliminary exploration of the use of the methods 
derived, to show the broad effects of experimental conditions, to evaluate the 
accuracy of the derived quantities, and to find the optimum conditions for 
accuracy in further work. 

As a result of a complete analysis of the transient curve, four parameters 
associated with thermal transfer from a hot solid to the surrounding liquid could 
be obtained. These include two thermal resistance terms, one associated with 
a purely conductive flow (r in the paper cited) and one associated with a shunt 
flow due to convection alone (R). The other two parameters derived are C, 
a capacitance term and probably a measure of the thermal capacity of the 
Stagnant film around the heated solid, and Z, an inductance term tentatively 


associated with the kinetic energy stored in the system of circulating fluid around 
the heated body. 


II. EXPERIMENTAL DETAILS 
The system selected for experimental study consisted of a heated resistance 
alloy wire which was held horizontally by stout current lead and immersed in 
the fluid under test. The wire was heated by an A.C. current from a source of 
steady e.m.f. through a variable resistance and current meter. The wires used 
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- varied in diameter from 0-610 to 0:274mm (23-32 S.W.G.) and the lengths 


varied from 1-5 to 2-5cm. The heat flux density for unit length of a wire 
heated by a current J of diameter d and specific (electrical) resistance p is given 
by 

SHON SHE batik F saagtteeran siesta Seely amine oe (1) 


with J the Joule’s equivalent. For a fixed experimental Set-up q is thus directly 
proportional to I?. 


The temperature difference between the wire and the fluid (6) was obtained 
by means of a thermocouple of wire at least 3 gauges finer than the heater wire. 
One junction of this couple was welded to the centre of the heater wire and the 
other (immersed in the fluid) was situated some 5mm away and in the same 
horizontal plane. The thermocouple was connected through a variable shunt 
box to a recording galvanometer. Relative measures of 6 can be obtained from 
the galvanometer reading and absolute readings can be obtained from calibration 
of a like couple in baths of different known temperatures. 


The liquids examined were water, aqueous solutions, and various organic 
liquids. Some 200ml was used in any test. No special purification was 
attempted. The vessel containing the liquid under test and the hot wire thermo- 
couple device were contained in an external thermostat at a temperature ranging 
from 12 to 70°C and the vessel was allowed to come to thermal equilibrium 
before any measure was attempted. A predetermined heating current was then 
switched on and the temperature difference recorder was allowed to run until the 
temperature became steady. When the time taken was in excess of 10 min an 
ordinary bench reflecting galvanometer with frequent reading could be used. 
Processes which could have been completed in less than 5 sec have not yet been 
examined, although there appears to be evidence that such fast transients could 
be obtained under certain conditions. Particularly quick acting recorders are 
necessary in the analysis of these transients and these have not yet been used in 
this field. 

The experimental factors which might be expected to affect the nature of 
the convective transient are the diameter of the wire, its temperature excess 0, 
and the thermal nature of the fluid, including the thermal conductivity and the 
convective modulus. Changes in 6 provide the easiest experimental variation, 
limitation being set only by the onset of boiling at the upper limit and the 
sensitivity of the galvanometer at the lower limit. However, as will be seen 
later, there is a great tendency for the time constants to become excessively 
long at low values of 6 and the ultimate limiting factors may be the stability of 
the heating current and of the thermostat bath. Changes in d rath the present 
form of device are limited to a comparatively narrow range. Below 32 8.W.G. 
it is difficult to get suitable thermocouples which do not produce an appreciable 
thermal disturbance at the centre of the heated wire. Wires greater than 
22 S.W.G. tend to draw an excessive heating current and seriously disturb 


external thermostating. The range of thermal conductivity in non-metallic 


‘liquids is comparatively restricted (Bosworth 1952). On the other hand there 


are readily available fluids covering a wide range of convective moduli; this 
ir ¥ 
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modulus is an experimental variable which may be readily changed by using 
viscous liquids at different compositions or temperatures. Mixtures used include 
the systems glycerol—water and sucrose—water. 


III. EVALUATION OF EXPERIMENTAL DATA 

As already intimated, relative values of the heat flux q and the temperature 
difference 8 may be obtained more readily than absolute figures. However, 
the first step in the deduction of the circuit parameters is concerned solely with 
dimensionless quantities, A and B, which are ratios of times and impedances 
respectively as derived from the transient curves. From A and B, using data 
from the tables and curves given in the previous paper (Bosworth and Groden 
1960), it is possible to derive the quantities X,Y, and R/r, all of which are 


TIME (MIN) 


Fig. 1—Transient curves, plotted as 6/J? against time. 


dimensionless. To proceed further it is necessary to take into account the 
absolute value of t, and then get the quantity +/(ZC), or alternatively to take into 
account the final thermal impedance Z, which is equal to Rr/(R-+-r) or 6/I? for 
any fixed experimental set up. From the known absolute (or relative) values of 
6/q it is then possible to derive R and r. From 


4 L IO == VR heed tern <8 uae (2) 


and from /(LC) it is then possible to evaluate both L and 1/C on the same (or 
relative) scale. 

As an illustration we will consider data derived from a wire of diameter 
0-05105 cm immersed in 60-4 per cent. w/w sucrose solution at 12°C. The 
heating current ranged from 0-230 A through 0-362, 0-540, 0-790, 1-220, to 
1-590 A. The transient curves, plotted in the form of 6/J? versus time, are 
shown in Figure 1. The lowest heating rate gave a curve with a very indistines 
maximum, at least within the time of the experiment recorded, while the curves 


t 
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associated with the higher heating rates showed transients passing successively 
faster and faster through a maximum and then approaching the final value from 
above. Some of the curves show signs of a severely damped oscillation after 
passing through the first maximum. 


From readings taken from the third curve we obtain the following four 
characteristics : 


De aca A=3-40 
t, 92 sec Bi res 
Li 20-4 arbitrary units ) Fe 

: oe ot B=1-415. 
VA 29-9 arbitrary units) 


From Figure 4 of the previous paper it is now possible to read off both 


XY=2-3 and R/r=1-80, 
from which 

Mg == 2 d (aay 

From the value of Z, (measured in arbitrary units as 6/7?) we then arrive 
at the values 
R=57-2 and r=31-7 arbitrary units. 
From the values of R and Y we obtain 
4/(L/C)=157 in the same arbitrary units, 


and finally using the value of t, and the ratio »/(LC)/t,, depending on X and 
reading off from Table 2 of the previous paper, we get 


4/(LC)=131 sec. 


Accordingly 
L=20600 arbitrary units, and 
C=0-83 reciprocal units. 


IV. EXPERIMENTAL RESULTS 
Analysis of the curves shown in Figure 1 may be carried out by the method 
outlined above. However, it is more convenient to convert the purely relative 
figures to absolute figures from the known diameter and electrical resistivity of 
the wire concerned, using equation (1). For a wire of diameter 0-05105 em 
with specific resistivity p=10-* Q cm, the heat flux per unit length becomes 


q=0-0711? cal em} sec, wee eee eee eee (3) 


so that in turn rR/(R-+r) becomes 
(6/q)anai=14 (0/2?) finai- 


The appropriate values of the circnit parameters derived from the six curves 
shown in Figure 1 are given in Table 1, while Table 2 contains corresponding 
- values from the same wire in different solutions and Table 3 for different wires 


_in various liquids. 
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TABLE 1 
THERMAL CIRCUIT PARAMETERS ASSOCIATED WITH A WIRE OF DIAMETER 0-0511 cm IN 60-4 PER CENT. 
SUCROSE SOLUTION AT 12 °C 


i) r R L 6 
(degC) (degC em sec cal-1) | (degC cm see cal-) | (degC cm sec~* eal-1)| (cal (degC)— em) 
1:7 540 2700 25,000,000 0-160 
3°15 465 1200 1,300,000 0-065 : 
5-9 445 800 280,000 0-060 I 
11-5 400 750 160,000 0-049 ; 
23°2 376 535 35,000 | 0-041 
34-4 364 410 12,000 | 0-054 
TABLE 2 
THERMAL CIRCUIT PARAMETERS ASSOCIATED WITH A WIRE OF 0-051] CM DIAMETER IN SUCROSE : 
SOLUTION 
Strength Bath 9 r | R L Cc ' 
of Temper- (degC) (degC (degC (degC (cal 
Solution ature cm sec cal—!) | cm sec cal-) | cm sec? cal-) | (degC)? cm=—") 
60:4% 16 °C 2:6 440 2300 1,700,000 0-072 
377 380 1700 860,000 0-049 
5:5 400 1260 410,000 0-045 
9-2 460 1200 170,000 0-044 — 
12:8 440 780 78,000 0-049 
18-5 410 700 53,000 0-036 
25°5 390 700 39,000 0-030 
34 390 490 11,000 0-045 
59 380 350 : 3,100 0-046 
63 440 280 2,500 0-035 
79 520 170 710 0-037 
86 620 63 510 0-048 
112 320 280 1,400 0-029 
138 245 190 1,100 0-046 
40-5°C | 3:4 650 2600 2,740,000 0-048 
7:6 460 590 49,000 0-043 
13-0 320 680 19,000 0-043 
42 290 360 2,900 0-041 
69 °C 5:6 420 1400 240,000 0-022 
W7-1 500 400 24,000 0-022 
24-6 430 300 2,000 0-037 
45-6 490 200 410 0-021 
81-2% 18 °C 132 520 520 12,800 0-050 
554% 19)°C 6-5 430 640 4,200 0-045 
40-0 340 250 220 0-054 


ratio recorded in the data used in Table 1, a 1p 
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TABLE 3 
THERMAL CIRCUIT PARAMETERS FROM VARIOUS SOURCES 


Wire | 
Diameter | Liquid fe: saline te , z ‘ 
q Per-| (degC) (degC em | (degC em | (degC em | (cal (degC)-1 
(mm) ature sec cal-) | sec cal!) | sec? cal-1) cm?) 
0-457 40% is) SG: 2:8 530 680 53,000 0-083 
sucrose 4:3 840 400 8,700 0-077 
solution* 
0-457 70% —- 47 440 1160 285,000 0-044 
sucrose* 
0-1016 water 22 °C 26 | 12-1 42-1 2,700 0-168 
0-1016 carbon — 4-0 90 330 280,000 0-084 
tetrachloride 
(?) 
0-1016 toluene —- 2-4 64 620 18,000 0-170 
Bio] 61 510 3,800 0-089 
4-4 65 350 3,500 0-092 
37 55 130 2,600 0-090 


* Data from published curves (Bosworth 1946a). 


V. ACCURACY 

Before considering possible relationships between the circuit parameters 
and the properties of the heated wire and the convective medium, some thought 
must be given to possible sources of error. Two obvious sources of inaccuracy 
immediately come to mind. The device used cannot strictly be described by 
any simple mathematical expression. The small thermocouple has been used to 
measure the temperature difference between the central portion of the hot wire 
and the surrounding fluid. In practice, of course, the presence of the couple 
alters the electrical conductivity of the central portion of the wire and so disturbs 
the rate of heat generation. Again, the arms of the couple carry heat away from 
the heated zone and so tend again to depress the temperature of the zone for 


which a temperature measurement is required. Disturbing influences may be 


reduced by using finer wires for the thermocouple and may become insignificant 
with sufficiently fine wire. However, in an initial broad study of the allied 
phenomena this particular aspect has not yet received quantitative attention. 


For any particular experimental set-up, assuming that the actual temperature 
is a representative central temperature of the wire, then the accuracy for all four 
parameters is largely dominated by the relationship between the A’s, or ratio of 
the times for attaining standard points on the transient, and the B’s, or ratio of 
the temperature differences. Figure 4 in the previous paper shows this relation 


ship. At a value of B equal to 1 -2, which is the average value of this experimental 
er cent. error in the ratio A will 
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give rise to 3 per cent. error in the ratio of R/r and thus subsequently in the 
derived parameters. Again, at a value of A equal to 2-5 (the average of the 
experimental figures used above) an error of 1 per cent. in the value of B will 
give rise to a 10 per cent. or larger error in the ratio of R/r and thus in subse- 
quently derived parameters. Accordingly, it is vitally important, if reliable 
circuit parameters are to be obtained, that the measurements of B should be made 
with the utmost care. It is indeed fortunate that the amplitudes of transient 
records can be recorded and measured with considerable accuracy—with a 
higher accuracy for example, than that which may be recorded for the time 
required to reach a turning point. In the latter case the flat nature of the 
maximum and the asymmetrical shape of the peak considerably limits the 
accuracy of assessment of the true turning point. However, assuming that we 
can measure Z, (or the amplitude at the peak) with a reasonably high accuracy 
the value of all subsequent data will primarily be limited by the accuracy with 
which we may measure Z, (or the final impedance). For this reason it is most 
desirable that current stabilizers should be used in the heating circuit. Any 
fluctuation in the rate of heating will vitiate attempts to attain a highly accurate 
measurement of B as Z,/Z,. Again, if the transient involves a particularly long 
time the experimentalist may read an apparent value of Z which has in faet not 
settled down to the final Z. Again, if the period is long and as a considerable 
source of internal heat must be applied to the fluid during the course of the 
transient, the apparent Z, may in fact not be the true Z, because the properties 
of the fluid have been changed. Long-period transients may therefore be 
expected to give false readings of the thermal circuit parameters. 


On the other hand, if the thermal transient is associated with too short a 
period the galvanometer or recorder may not follow truly the whole course of the 
transient, particularly as the low resistance associated with the thermocouple 
circuit may cause serious overdamping of the galvanometer. Errors due to 
this cause usually first show up in abnormally high values of ¢, (the time taken 
first to cross Zp) relative to the value of t,. Overdamping of the galvanometer 
gives a spuriously low value of A, which has shown up in many instances in 
which a bench galvanometer has been replaced by a quick-acting potentiometric 
recorder which gives an apparently much higher value of A for the same conditions 
as those studied with a reflecting galvanometer. Generally speaking therefore 
the method described in these two papers can be applied to very quick-acting 
transients only if suitable quick-acting recorders are used. 


A study of the effect of errors due to mis-estimation of the final value of ae 
may be carried out on curve A of Figure 1. The final value of Z obviously 
is attained outside the range of time recorded in the figure given (that is, up to 
10 min) and has to be taken as 32 arbitrary units. Had we postulated that at the 
time of 30 min Z had in fact not yet reached its final value and that instead a 
longer time would have given Z,=31, then the values of Y and Y which had been 
computed as 11-6 and 4-65 would instead be 14-2 and 10 -3, with the associated 
circuit parameters in the order given in Table 1: 580, 1620, 42 «108, and 0-150. 
Conversely, an assumption that due to heating of the fluid the value of 32 had 


the wire is 
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already overshot the final true value and Z,=33 is a better representation, then 


the four parameters would have been estimated as at 510, 5050, 14108, and 
0-185. 

The apparent abnormally high values of 0, R, and L at the lowest heating 
rate cannot therefore be accepted at this stage and a more carefully designed 
experiment appears to be required. 


VI. CONCLUSIONS FROM EXPERIMENTAL DATA 

The variation in the four circuit parameters, as shown in Tables 1-3, while 
subject to fairly large errors, particularly at the lowest heating rate, are sufficiently 
pronounced to make a series of broad generalizations. 

All the four properties appear to change somewhat with change in the 
nature of the convective liquid and with the size of the heating surface. However, 
the most pronounced change is that associated with the change in the overall 
temperature difference 0. The two properties r and CO, those which have 
tentatively been associated with the conductive contribution to the total heat 
flux, are practically unchanged with change of 0. On the other hand, the 
properties & and Z, tentatively attributed to the convective contribution to the 
total heat flux, both decrease rapidly with increase of 0, the change in L being 
particularly pronounced. This behaviour, of course, is, at least in its general 
characteristics, in accordance with the well-known features of convective transfer. 
However, a quantitative examination of the magnitude of the quantities con- 
cerned could add considerable knowledge to this subject. 

The capacity term, for example, appears to vary with the diameter of the 
wire concerned. An average figure of 0-050 cal em—! (degC)-! was obtained 
for a wire of 0-511 mm diameter. This figure may be compared with the thermal 
capacity per centimetre of the wire itself. The corresponding capacity C7? of 


01= 1npcd? cal em—} (degC)-1, 


where 0 is the density of the wire, ¢ its specific heat, and d its diameter. 


Inserting numerical values we get 
01=1-7 x10-3 cal em—! (degC)-}, 


which is only about 1/30 of the measured value (C) of the capacity per unit 
length of the wire. The thermal capacity of whatever is contributing to the 
effective GC as deduced from the nature of the thermal transient is thus much 
larger than the thermal capacity of the heated solid alone. The thickness of 
the fluid layer surrounding the wire required to give a capacity of 0-05 
cal cm-! (degC)-1, using the Janovsky and Archangelsky (1928) data for the 
volumetric specific heat of sucrose solution, amounts to 1-1 mm. Taking the 
capacity of the toluene films surrounding a wire of 0-1016 mm diameter as 


0-09 cal em! (degO)—! (Table 3) gives a corresponding sheet thickness of 1-0 mm. 
- These figures may be compared with the estimated film thickness of 1 mm in 
‘mild thermal convection given by Hckert and Soehngen (1952). 
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The value of the parameter r in Tables 1-3 must be associated with the flow 
of heat from the wire before a condition of convection sets in. This involves 
an unsteady state problem in heat flow in a cylindrical system and is best not 
exhaustively treated until more complete and more accurate data have been 
compiled. Suffice it to say, however, that the solution will depend only on the 
diameter of the wire (possibly with correction due to the presence of the thermo- 
couple) and on the thermal conductivity of the liquids. Analysis thus offers a 
possibility of eliminating the disturbing effects of convection on the measurements 
of conduction in mobile liquids, and this technique may well lead to a quick 
laboratory method for comparing the thermal conductivities of mobile liquids, 

The factors R, which are highly sensitive to the temperature excess 6, are 
in turn indices to the convective moduli. As such these factors describe the 
ways in which the additional heat flow rate varies with the overall temperature 
difference in natural thermal convection. 

The factors D were, in the first instance, thought to be interpreted as the 
measure of the kinetic energy stored in the system of circulating convection 
currents per unit of squared thermal flux. One is tempted to write, on analogy 
with the corresponding expression in electrodynamics, for the kinetic energy K, 


Kes lpe 2 een (4) 


This, however, on the basis of the definitions and units so far used, is dimensionally 
incorrect in just the same way as would the expression 


V=4( | aty?/C 


be dimensionally incorrect for the thermal (potential) energy stored in the 
capacity. 

These formal difficulties could, according to the author (Bosworth 19465), 
be neatly avoided by a redefinition of the thermal flux as a definition of a flow of 
entropy. Equation (3) then becomes 


0-071 
AS 


J3 BU oom =) 865k") Gays cele ees (5) 


with 7 the absolute temperature of the sink. On the basis of the units previously 
used, namely, joules (degC)-! em—}, q becomes 


0-30 
qI= al * thermal amperes/em. .......... (6) 


On this basis the various properties listed in Tables 1-3 could then be changed 
to new units as thermal ohms.cm (or (deg)? sec em joules!) for the second 
and third columns, thermal henry em (or (degC)? sec? em joules—!) for the fourth 
column, and thermal farads per cm (or joules em-! (degC)-*) for the last. The 


second, third, and fourth columns should be multiplied by T/4-2 and the fifth 
by 4:2/T7. 
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Equation (4) now becomes dimensionally correct and K is measured in 
joules per cm of the heated wire. As we are concerned with the convective 
contribution to the heat current we put 


q=6/R 
or 
K=3L6?/R?. 
To find the way in which the kinetic energy changes with the temperature 


difference 6, it is therefore necessary to find the relationship between L/R? and 0. 
Data from the Tables 1 and 2 have been plotted in Figure 2 in the form log L/R? 


Loe 0) 


3-5 To 1-5 


LOG tf/R2 
Fig. 2.—Relation between L/R* and 90. 


-yersus log 6. The differently marked points corresponding to different bath 
temperatures fall on different parallel straight lines, all having the slope —1-0. 


We accordingly put 
L/R?=a0-1, 


with the factor a as a constant for any particular wire in any particular fluid. 
se: * 2 ee ee ees ey (8) 
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so that it appears that the kinetic energy of the convective current increases 
linearly with the temperature excess. For the 0-511 mm wire in 60-4 per cent. 
sucrose solution at 12 °C the quantity a takes the value 0-041 J (degC)-* em~?, 
and at 16°C, 0-025 J (degC)! cm-}!. Accordingly, A=0-020 J/em at 12 a BE 
and 0-012 J/em at 16 °C, for a temperature rise of 1 degC and with somewhat 
lower figures in the less viscous fluids. Now the total volume of fluid in which 
the hot wire was immersed was some 200 ml, and at the most some 30 ml/em 
of wire could be involved in the system of convection currents. The expected 
kinetic energy in such a system could only be of the order of ergs/em. The 
phenomenon giving rise to the inductive features of thermal convection is of 
one or two orders of magnitude larger than that expected from the kinetic energy 
of the system of convection chimneys. 


Presumably any process in which the heat flow rate is partially conditioned 
by states of the system arising from absorption of energy in a non-thermal form 
from the initial heat flow will add an inductive-like element to the thermal 
circuit. The storage of heat in a process of surface evaporation would act in 
this way. In some instances cited in Table 2 the liquid was actually caused to 
boil on the wire (although the liquid medium as a whole did not boil). While 
in these instances the “ inductance ” was larger than the extrapolated value from 
the non-boiling region, however, the change was not one of a different order of 
magnitude, as might have been expected. Accordingly, there arises the suggestion 
that the unexpectedly large values associated with thermal induction without 
surface boiling may be due to surface chemical phase change, for example, to 
deaeration of the liquid near the wire and subsequent re-solution of microscopic 
bubbles in the slightly cooler region. This suggests that thermal transients 
could profitably be studied after complete deaeration of fluids under test. 


As an alternative interpretation of the high value of thermal inductance 
there may be cited our general lack of complete knowledge of the mechanism 
involved in the passage of layers of different fluid density over a solid surface— 
the so-called tea-pot effect (Reiner 1956). The tendency of flowing streams of 
different density to adhere to solid surfaces appears to be several orders of 
magnitude greater than that expected from normal surface (tension) phenomena. 
Reiner has suggested that the creation of micro-eddies is responsible. If upheld, 
a Similar mechanism may also be the cause of the apparently high kinetic energy 
term connected with the act of setting up the convection chimney. 
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SHORT COMMUNICATIONS 


AN ANALOGY BETWEEN THE EFFECTS OF FLUCTUATING ELECTRIC 
FIELDS AND STEADY MAGNETIC FIELDS IN ISOTROPIC CONDUCTORS 
WHEN A UNIVERSAL RELAXATION TIME CANNOT BE DEFINED* 


By E. J. Mooreyt 


It is well known that when a “ universal ” time of relaxation (zt) exists, the 
influence of a harmonically varying electric field (Foe!) on the transport 
properties of a solid may be taken into account by replacing t by 7t/(1+iwr). 
Dingle (1956a) demonstrated that, for an isotropic solid, the effect of a steady 
magnetic field may similarly be obtained by replacing t by 7/(1+jQc) with an 
applied d.c. electric field, and by t/[1+(iw+jQ)r] with an a.c. field. (Here 
j?=—1, ij A—1, and Q=(—e)H/me is the circular frequency of precession of an 
electron.) The object of the present note is to show that this analogy between 
a high frequency electric field and a steady magnetic field still exists, even when a 
“ universal ’’ relaxation time cannot be defined. 

For ease of exposition, we consider a one-band model and assume that the 
applied fields are an electric field F=(P,,F,,0) cel’ and a constant magnetic 
field H=(0,0,H). The acceleration of a charge carrier (charge —e, effective 
mass m, energy H=4(h/m)? | k |?, wave vector k, and velocity v=hk/m) is then 


: € ih 
As usual, we start with the Boltzmann equation 
Of if ) 2 
eae “INAS | JON | nl Reng nec. O.9 om Sb. 1 
ot <(F+ oe / Ve Ot J cont.” Le 


and assume a first-order solution of the form (Wilson 1953, p. 210) 
Ofo 
Se ede CUE ee oss cow habe c cists 2 


Here, f is the electron distribution function, [0f/0t].on. is the rate of change of f due 
to the various scattering mechanisms, f, is the equilibrium distribution function, 
and c(#) is a function of the electron energy and the applied fields. 

| Now, for the more important of those scattering mechanisms for which a 
universal relaxation time does not exist—acoustic phonon scattering in metals 
below the Debye temperature (Wilson 1953, p. 263) and optical phonon scattering 
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in polar semiconductors (Howarth and Sondheimer 1953)—the collision term 
may be written 

| =k-L(c) so taal on ra ees eae (3) 

at coll. 
where c= (¢,¢,0), L(c)=(L[e,],L[¢,],0), and L is a linear operator, i.e. if « and 6 
are any functions independent of the electron energy, then 

Liac+8d)=abi(c)--Bii(d). .. viene one (4) 


Substituting (2) and (3) in (1) we find, correct to linear terms in the electric 
field, that the components of c satisfy the equations 


OC, ch 


L(¢,)+ at —Oc,= — me® Rene ee? (5) 
0c ch 
L(e,)-+ a +Qe¢,=— ae vie take babe (6) 
Now, by introducing a second imaginary j (j?=—1, ij 4—1) equations (5) and 
(6) may be combined into the single equation 
Coe eh 
L(C)+ = +jQ¢=— —-F, aug tS eee (7) 


where C=c,+jc,, and F=F,+jF,. 


It is now convenient to introduce a formal relaxation time t(@,Q) defined by 
ch 
C= — -Fr(o, 2). Sis, ona aie atte ee (8) 


Since F oce’o!, equation (8) is a solution of (7) provided +(w,Q) satisfies 


L[t(,Q)] +(i@ +jQ)t(w,Q)=1. .......... (9) 
With o=Q=0, equation (9) reduces to 
Li[r(6,0) |=) ee (10) 


where 7(0,0) corresponds to one of the formal relaxation times introduced by 
Dingle (1956b). 

Solution of the Boltzmann equation (9) is extremely difficult. The object 
of the present note is to direct attention to one interesting simplifying feature— 
that the solutions of (9) for the three field combinations (#=0, 20), 
(o40, Q=0), and (w40, QO) are essentially the same, so that soluticn of 
the Boltzmann equation for any one of the above cases in fact implies a full 
solution of the problem for all the cases. 

Finally, it is of interest to compare the “ universal’? and the “ formal ” 
relaxation time theories. As mentioned above, the influence of a uniform 
magnetic field and a high frequency electric field on the transport properties 
may be obtained by replacing t—! by t-1+iw+jQ. Analogously, for the formal 
relaxation time there is the replacement of the collision operator L by L+io+jQ. 


SHORT COMMUNICATIONS 97 


The author is indebted to the Commonwealth Scientific and Industrial 
Research Organization for the award of a Studentship. 


References 

Dinetz, R. B. (1956a).—Physica 22: 701. 

Dinetz, R. B. (19566).—Physica 22: 698. 

Howarru, D. J., and Sonpurrmer, E. H. (1953).—Proc. Roy. Soc. A 219: 53. 
Witson, A. H. (1953).—‘‘ The Theory of Metals.” 2nd Ed. (Cambridge Univ. Press.) 


AUSTRALASIAN MEDICAL PUBLISHING CO. LTD. 
SEAMER AND ARUNDEL STS., GLEBE, SYDNEY 


NOTICE TO CONTRIBUTORS 


‘ GENERAL.—Papers will be considered for publication irrespective of the 
organization to which the authors are attached. Intending contributors, 

in preparing papers for submission, should follow the general style adopted 

a ise os oe A “Quide to Authors”? may be obtained on application 
o the or. 


Papers, which should be written as concisely as possible, should be double- 
space typed on one side of the paper only and with liberal side margins. 
Every page, including those carrying tables and illustrations, should be 
numbered. _ Lengthy manuscripts should be accompanied by a table of 
contents with headings appropriately numbered. The organization to which 
the author is attached should be indicated in a footnote. 


The original typescript and one carbon copy should be forwarded. 


- Summary.—A short summary, which may be paragraphed, should precede 
the introduction to the paper. It should be of an informative character and 
written in such a way that it can be used by abstracting journals without 
amendment. 


. REFERENCES.—References are to be cited in the text by the year of 
publication, e.g. Thomson (1948), and are to be arranged alphabetically, 
giving the author’s name and initials followed by the year of publication, 
title of the paper (if desired), title of the periodical, volume, and pages, thus : 


Thomson, G. P. (1948).—Proc. Phys. Soc. Lond. 61: 403-16. 


Abbreviations of titles of periodicals should conform to those used in “ A 
World List of Scientific Periodicals ”’. 


. ILLUSTRATIONS.—Line drawings should be made with indian ink on white 
drawing paper (preferably bristol board), tracing paper, or graph paper. 
Only blue or very pale grey ruled graph paper should be used as yellow, green, 
and red lines are difficult to screen out. Blue linen should not be used. Lines 
should be sufficiently thick and symbols sufficiently large to allow of reduction 
to a page size of 73 by 5 in. (or less). It is suggested that figures be drawn 
approximately twice the width at which they are to be printed. Lettering 
and numbering should only be pencilled on the drawings as the appropriate 
type will be substituted before the blocks are made. Original drawings 
should be submitted but, where possible, photographic copies or blue-prints 
of the originals should also be forwarded. 


Half-tone photographs, which should be included only where essential, should 
be on glossy paper and show as much contrast as possible. 


. REPRINTS.—Individual authors receive 40 copies of each paper free, two 
authors receive 30 each, and three or more authors receive 20 each, but 


additional copies may be purchased. 


CONTENTS 


PAGE 


General Method of Exact Solution of the Concentration-dependent Diffusion 
Equation. By J. R. Philip 


The Function inverfe 0. By J. R. Philip 


The Distribution of Ions formed by Attachment of Electrons moving in a 
Steady State of Motion through a Gas. By C. A. Hurst and L. G. H. 
Huxley a a7 ie i <a ee Pe a 


The Design of Photographic Objectives of the Triplet Family. II. The 
Initial Design of Compound Triplet Systems. By F. C. Cruickshank 


The Damping of Water Waves by Surface Films. By R. G. Vines 


Sunrise and Eclipse Effects on the Ionosphere at Brisbane. By G. G. 
Bowman lt - ~ ; ae ee 


A Relationship between Spread-F and the Height of the F, Ionospheric 
Layer. By G. G. Bowman = =e : me aa 


Thermal Transients associated with Natural Convection. By R. ©. L. 
Bosworth and C. M. Groden 


Attempts to measure the Inductive Element associated with the Natural 
Convection of Heat. By BR. C. L. Bosworth .. 


Short Communications 


An Analogy between the Effects of Fluctuating Electric Fields and Steady 
Magnetic Fields in Isotropic Conductors when a Universal Relaxation 
Time cannot be Defined. By E. J. Moore 


-_—ooOo 
AUSTRALASIAN MEDICAL PUBLISHING CO. LTD. 
SEAMER AND ARUNDEL STS,, GLEBE, SYDNEY 


13 


69 


73 


84 


95 


—— — 


